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ABSTRACT 



o 

5_( . Context. The Cepheid period-luminosity (PL) relation is unquestionably one of the most powerful tools at our disposal for determin- 

"^ ' ing the extragalactic distance scale. While significant progress has been made in the past few years towards its understanding and 

rrt ' characterization both on the observational and theoretical sides, the debate on the influence that chemical composition may have on 

the PL relation is still unsettled. 

Aims. With the aim to assess the influence of the stellar iron content on the PL relation in the V and K bands, we have related the 
04 ' V-band and the A'-band residuals from the standard PL relations of Freedman et al. (2001) and Persson et al. (2004), respectively, to 

> ; [Fe/H]. 

^^ , Methods. We used direct measurements of the iron abundances of 68 Galactic and Magellanic Cepheids from FEROS and UVES 

0^ ■ high-resolution and high signal-to-noise spectra. 

Results. We find a mean iron abundance ([Fe/H]) about solar (cr = 0.10) for our Galactic sample (32 stars), — 0.33 dex (cr = 0.13) 

for the Large Magellanic Cloud (LMC) sample (22 stars) and ~ -0.75 dex {cr = 0.08) for the Small Magellanic Cloud (SMC) sample 

(14 stars). Our abundance measurements of the Magellanic Cepheids double the number of stars studied up to now at high resolution. 

The metallicity affects the K-band Cepheid PL relation and metal-rich Cepheids appear to be systematically fainter than metal-poor 

ones. These findings depend neither on the adopted distance scale for Galactic Cepheids nor on the adopted LMC distance modulus. 

Current data do not allow us to reach a firm conclusion concerning the metallicity dependence of the A'-band PL relation. The new 

Galactic distances indicate a small effect, whereas the old ones support a marginal effect. 

Conclusions. Recent robust estimates of the LMC distance and current results indicate that the Cepheid PL relation is not Universal. 
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1 . Introduction =0.02), independent of wavelength. Nonlinear convective mod- 

, , , , , ^ , . , T, . , els (e.g. Bono et al. 1999; Caputo et al. 2000; Caputo 2008) in- 

Smce the dawn of modern astronomy the Cepheid Period- ^^^^^ ^.^^.^^ ^ j dependence on the same interval of metal- 

Luminosity (PL) relation is a key tool in determining Galactic j.^. ^j^^ ^^ -^ q 4 -^ y^ q 3 .^ j ^^^ ^^ mag in 

and extragalactic distances In spite of its fundamental impor- ^^ ^^^ ^^ j^ ^p^ ^ ^ Moreover, the predicted change in these 

tance, the debate on the role played by the chemical compost- j^^^^^ ^^^^j^ -^ ^^^^ ^^^^ metal-rich Cepheids are fainter than 

tion on the pulsation pi-operties of Cepheids is sti 1 opeii, with ^^^^^1 , ^nes, at variance with the results of the linear models, 

different theoretical models and observational results leading to pi^rentino et al. (2002) and, more recendy, Marconi, Musella & 

markedly different conclusions. ,,.,.. Fiorentino (2005) investigations, also based on nonlinear mod- 

From the theoretical point of view pulsatioii models by dif- ^j^^ ^^^^ ^^^^.^ ^^ ^1^^ ^ dependence on the helium 

ferent groups lead to substantially different results. Linear mod- ahnnHancp 
els (e.g. Chiosi et al. 1992; Sandage et al. 1999; Baraffe & 

Alibert 2001), based on nonadiabatic radiative models, suggest a ^n the observational side, the majority of the constraints 

moderate dependence of the PL relation on the metallicity. The '^""^^s ^om indirect measurements of the metallicity, mosdy in 

predicted change at log(P) = 1 is less than 0.1 mag for metal ^^ternal galaxies, such as oxygen nebular abundances derived 

abundances ranging from the SMC (Z=0.004) to the Galaxy (Z ^om spectra of H II regions at the same Galactocentric distance 

as the Cepheid fields (e.g. Sasselov et al. 1997; Kenmcutt et al. 

Send ojfprint requests to: G. Bono, e-mail: bono@mporzio . astro . it 1998; Sakai et al. 2004). These analyses indicate that metal-rich 

* Based on observations made with ESO Telescopes at Paranal and Cepheids are brighter than metal-poor ones (hence at variance 

La Silla Observatories under proposal ID 66.D-0571 with the predictions of nonlinear convective models), but it is 
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important to note that the results span a disappointingly large 
range of values (see Table 1 and Fig. 1). 

More recently Macri et al. (2006) found, by adopting a large 
sample of Cepheids in two different fields of NGC 4258 and the 
[O/H] gradient based on H II regions provided by Zaritsky et 
al. (1994), a metalhcity effect of y = -0.29 + 0.09 mag/dex. 
This galaxy has been adopted as a benchmark for estimating the 
metallicity effect, since an accurate geometrical distance based 
on water maser emission is available (Herrnstein et al. 2005). 
However, in a thorough investigation Tammann et al. (2007) 
suggested that the flat slope of the Period-Color relation of the 
Cepheids located in the inner metal-rich field could be due to 
a second parameter, likely helium, other than the metal abun- 
dance. Furthermore, Bono et al. (2008) found, using the new and 
more accurate nebular oxygen abundances for a good sample of 
H II region in NGC 4258 provided by Diaz et al. (2000), a shal- 
lower abundance gradient. In particular, the new estimates seem 
to suggest that both the inner and the outer field might have a 
mean oxygen abundance very similar to LMC. They also found 
a very good agreement between predicted and observed Period- 
Wesenheit {V, I) relation. Nonlinear convective models predict 
for this relation a metallicity effect of y = +0.05 ±0.03 mag/dex. 

On the basis of independent distances for 18 galaxies based 
on Cepheid and on the Tip of the Red Giant Branch, Tammann 
et al. (2007) found a smafl metaflicity eff^ect (y = -0.017 + 0.1 13 
mag/dex). On the other hand, Fouque et al. (2007) using a sam- 
ple of 59 Galactic Cepheids whose distances were estimated us- 
ing different methods - HST trigonometric parallaxes (Benedict 
et al. 2007), revised Hipparcos parallaxes (van Leeuwen et al. 
2007), infrared surface brightness method (Fouque & Gieren 
1997), and interferometric Baade-Wesselink method (Kervella 
et al. 2004), zero-age-main-sequence fitting of open clusters 
(Turner & Burke 2002) - found no significant difference be- 
tween optical and Near-Infrared (NIR) slopes of Galactic and 
LMC Cepheids (Udalski et al. 1999; Persson et al. 2004). 

An alternative approach is to measure directly the metal con- 
tent of Cepheid stars, which, so far, has been attempted only 
by few studies, primarily focused on stars of our own Galaxy 
(Luck & Lambert, 1992; Fry & Carney 1997; Luck et al. 1998; 
Andrievsky et al. 2002a,b,c; Luck et al. 2003; Andrievsky et 
al. 2004). Fry & Carney (1997, hereafter FC97), for instance, 
have derived iron and a-element abundances for 23 Galactic 
Cepheids from high resolution and high signal-to-noise spectra. 
They found a spread in [Fe/H] of about 0.4 dex, which they claim 
is real. Using approximately half of their sample, the stars be- 
longing to clusters or associations, they have made a preliminary 
evaluation of metallicity effects on the zero point of the PL re- 
lation, finding that metal-rich Cepheids are brighter than metal- 
poor ones. Thus, finding a result similar to the studies based on 
indirect measurements of the metallicity. 

The impressive observational effort carried out by 
Andrievsky and collaborators (Andrievsky et al. 2002a, 2002b, 
2002c; Luck et al. 2003; Andrievsky et al. 2004; Kovtyukhet al. 
2005b) has, instead, taken advantage of high resolution spectra 
of 130 Galactic Cepheids (collected with different instruments 
at different telescopes) in order to determine their chemical 
composition and study the Galactic abundance gradient. The 
sample covers a range of Galactocentric distances from 4 to 14 
kpc. The emerging picture can be best described by a relatively 
steep gradient (about -0.14 dex kpc"') for Galactocentric 
distances less than 7 kpc, followed by a much shallower slope 
(~ -0.03 dex kpc"') between 7 and 10 kpc, a discontinuity at 
approximately 10 kpc and a nearly constant metallicity of about 
-0.2 dex towards larger Galactocentric distances, out to about 
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Fig. 1. Comparison of recent results for the metallicity sen- 
sitivity of Cepheid distances. FM1990: Freedman & Madore 
(1990); G1994: Gould (1994); Kol997: Kochanek (1997) 
S1997: Sasselov et al (1997); Kel998: Kennicutt et al (1998) 
F2001: Freedman et al (2001); U2001: Udalski et al (2001) 
C2002: Ciardullo et al (2002); Sa2004: Sakai et al (2004) 
Sto2004: Storm et al (2004); Gro2004: Groenewegen et al 
(2004); M2006: Macri et al. (2006); Sah2006: Saha et al. (2006); 
Be2007: Benedict et al. (2007); Fo2007: Fouque et al. (2007); 
Tam2007: Tammann et al. (2007); B2008: Bono et al. (2008). 
See Table 1 . 



14 kpc. In relation to our work, it is important to note that 
Andrievsky and collaborators did not investigate the effects of 
the chemical composition on the Cepheid PL relation: on the 
contrary, they used the PL relation to determine the distances of 
their stars. 

Outside our Galaxy, Luck & Lambert ( 1 992, hereafter LL92) 
have studied 10 Cepheids in the Magellanic Clouds (MCs). Five 
are in the Large MC (LMC) and five in the Small MC (SMC). 
For the former sample, they found a mean [Fe/H] of -0.36 dex 
with a dispersion of 0.3 dex, while for the latter one the mean 
[Fe/H] is -0.60 dex with a rather small dispersion of less than 
0.15 dex. A more recent study by Luck et al. (1998, hereafter 
L98) on 10 LMC Cepheids and 6 SMC Cepheids, 4 of which 
in common with LL92, confirmed the mean [Fe/H] value in 
the LMC (-0.30 dex), found very little evidence of a significant 
metallicity dispersion in the LMC (contrary to LL92, but sim- 
ilarly to the SMC), and slightly revised downwards the mean 
[Fe/H] of the SMC (-0.74 vs -0.60 found by LL92). 

Finally, we mention two studies that followed slightly differ- 
ent approaches. Groenewegen et al. (2004) have selected from 
the literature a sample of 37 Galactic, 10 LMC and 6 SMC 
Cepheids for which individual metallicity estimates and BVIK 
photometry were known. Their work aimed at investigating the 
metallicity dependence of the PL relation using individual metal- 
licity determinations as well as good individual distance esti- 
mates for Galactic Cepheids. They inferred a metallicity effect 
of about -0.27 ± 0.08 mag/dex in the zero point in VIWK, in the 



Romaniello et al.: The impact of chemical composition on Cepheid properties 3 

Table 1. Overview of recent results for the metallicity sensitivity of Cepheid distances. In the first column is listed the variation 
of the distance modulus jj. per dex of metallicity, the negative sign indicates that the true distance is longer than the one obtained 
neglecting the effect of the metallicity. In the second column is listed the elemental abundance used as reference for the metallicity. 
The third and fourth columns give the method and the reference of the different studies. See also Fig. 1. 
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[O/H] 
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[O/H] 
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[O/H] 





[Fe/H] 





[O/H] 


-0.24 ± 0.05 


[O/H] 


-0.21 ±0.19 


[Fe/H] 


-0.23 + 0.19 


[Fe/H] 


-0.29 + 0.19 


[Fe/H] 


-0.27 + 0.08 


[Fe/H] 


-0.39 + 0.03 


[Fe/H] 


-0.29 + 0.09 


[O/H] 


-0.10 + 0.03 


[Fe/H] 


-0.017 + 0.113 


[O/H] 





[Fe/H] 


+0.05 + 0.03 


[Fe/H] 



Analysis of Cepheids in 3 fields of M31 (BVRI bands) 

Comparison of Cepheids from 3 fields of M3 1 and LMC (BVRI bands) 

Simultaneous solution for distances to 17 galaxies (UBVRIJHK hands) 

Comparison of EROS observations of SMC and LMC Cepheids (VR 

bands) 

Comparison of HST observations of inner and outer fields of MlOl 

Comparison of 10 Cepheid galaxies with Tip of the Red Giant Branch 

distances 

Value adopted for the HST Key Project final result 

OGLE result comparing Cepheids in IC 1613 and MC (VI bands) 

Comparison of Planetary Nebula luminosity function distance scale and 

Surface Brightness fluctuation distance scale 

Comparison of 17 Cepheid galaxies with Tip of the Red Giant Branch 

distances 

Baade-Wesselink analysis of Galactic and SMC Cepheids (VK bands) 

Baade-Wesselink analysis of Galactic and SMC Cepheids (/ band) 

Baade-Wesselink analysis of Galactic and SMC Cepheids (W index) 

Compilation from the literature of distances and metallicities of 53 

Galactic and MC Cepheids (VIWK bands) 

Cepheid distances to SNe la host galaxies 

Cepheids in NGC 4258 and [O/H] gradient from Zaritsky et al. (1994) 

Weighted mean of Kennicutt, Macri and Groenewegen estimates 

Comparison between Cepheid and TRGB distances for 18 galaxies 

Comparison between the slopes of Galactic and LMC Cepheids 

Predicted Period- Wesenheit (V, I) relation 



Freedman & Madore (1990) 
Gould (1994) 
Kochanek (1997) 
Sasselovetal. (1997) 

Kennicutt et al. (1998) 
Kennicutt etal. (1998) 

Freedman et al. (2001) 
Udalski etal. (2001) 
CiarduUo et al. (2002) 

Sakai et al. (2004) 

Storm et al. (2004) 
Storm et al. (2004) 
Storm et al. (2004) 
Groenewegen et al. (2004) 

Saha et al. (2006) 
Macri et al. (2006) 
Benedict et al. (2007) 
Tammann et al. (2007) 
Fouque et al. (2007) 
Bono et al. (2008) 



sense that metal-rich Cepheids are brighter than the metal-poor 
ones (see Table 1 and Fig. 1, for a comparison with other stud- 
ies). Also Storm et al. (2004) discussed the effect of the metal- 
licity on the PL relation using 34 Galactic and 5 SMC Cepheids, 
for which they determined accurate individual distances with the 
Baade-Wesselink method. Assuming an average abundance for 
the SMC Cepheids of [Fe/H]= -0.7 and solar metallicity for the 
Galactic ones, they determined, in a purely differential way, the 
following corrections: -0.21 + 0.19 for the V and K bands, -0.23 
+ 0.19 for the /-band and -0.29 + 0.19 for the Wesenheit index 
W. These agree well with Groenewegen et al. (2004). 

Despite these ongoing observational efforts, it is important 
to underline that none of the observational studies undertaken 
so far has directly determined elemental abundances of a large 
sample of Cepheids in order to explicitly infer the metallicity 
effect on the PL relation, taking advantage of a sample that has 
been homogeneously analysed. 

The novelty of our approach consists exactly in this, i.e. in 
the homogeneous analysis of a large sample of stars (68) in three 
galaxies (the Milky Way, the LMC and the SMC) spanning a 
factor of ten in metallicity and for which distances and BVJK 
photometry are available. Preliminary results based on a sub- 
sample of the data discussed here were presented by Romaniello 
et al. (2005). Here we present the results about the iron content 
for the complete sample. In a forthcoming paper we will discuss 
the a-elements abundances. 

The paper is organised as follows. The data sample is pre- 
sented in Section 2. In Section 3 we thoroughly describe the data 
analysis and how we determine the metallicity of our stars. We 
compare our iron abundances with previous results in Section 
4. The dependence of the PL relation on [Fe/H] is discussed 



in Section 5. Finally, Section 6 summarizes our concluding re- 
marks. 



2. The data sample 

We have observed a total of 68 Galactic and Magellanic 
Cepheid stars. The spectra of the 32 Galactic stars were col- 
lected at the ESO 1.5m telescope on Cerro La Silla with 
the Fibrefed Extended Range Optical Spectrograph (FEROS, 
Pritchard 2004). Two fibres simultaneously feed the spectro- 
graph: one fibre records the object spectrum, while the second 
one is fed either by the sky background or by a calibration lamp 
to monitor the instrument stability during the exposure (the sky 
background was chosen in our case). The CCD is a thinned, back 
illuminated detector with 15 micron pixels (2048 x 4096 pixels). 
The resolving power is 48,000 and the accessible wavelength 
range is from 3,700 to 9,200 A. 

The spectra of the 22 Cepheids in the LMC and the 
14 Cepheids in the SMC were obtained at the VLT-Kueyen 
telescope on Cerro Paranal with the UV- Visual Echelle 
Spectrograph (UVES, Dekker et al. 2000; Kaufer et al. 2004) 
in service mode. The detector in the Red Arm is a mosaic of two 
CCDs (EEV + MIT/LL) with 15 micron pixels (2048 x 4096 
pixels). The resolving power is about 30,000 (corresponding to 
a slit of 100) and the spectral range covered in our spectra is 
from 4800 to 6800 A (580 nm setting). 

The 2-D raw spectra were run through the respective instru- 
ment Data Reduction softwares, yielding 1-D extracted, wave- 
length calibrated and rectified spectra. The normalization of the 
continuum was refined with the IRAF task continuum. The 1-D 
spectra were corrected for heliocentric velocity using the rvcorr 
and dopcor IRAF tasks. The latter was also used to apply the ra- 
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Table 2. Pulsation phases (0) and intrinsic parameters of the Galactic Cepheids. Both AP Pup and AX Vel were not included in the 
analysis of the metallicity effect because accurate distance estimates are not available in literature. In particular, for AP Pup we only 
listed the apparent mean magnitudes. In the last column is listed the duplicity status according to Szabados (2003): B - spectroscopic 
binary, Be - spectroscopic binary that needs confirmation, O - spectroscopic binary with known orbit, V - visual binary 



ID 


\0gP 


4> 


l^Old 


E(B - V)ou 


f^New 


E(B - \/)a,„„ 


Mb 


Mv 


Mk 


Duplicity 


ICar 


1.5509 


0.580 


8.99" 


0.170" 


8.56" 


0.147" 


-4.17 


-5.28 


-7.53 




UCar 


1.5891 


0.490 


10.97"' 


0.283"' 


10.87" 


0.265" 


-4.50 


-5.41 


-7.44 


B 


VCar 


0.8259 


0.375 


9.84'- 


0.174* 


10.09" 


0.169" 


-2.54 


-3.24 


-4.86 


B 


WZCar 


1.3620 


0.745 


12.92"' 


0.384"' 


12.69" 


0.370" 


-3.80 


-4.58 


-6.52 




VCen 


0.7399 


0.155 


9.18'' 


0.289" 


8.91" 


0.292" 


-2.41 


-2.99 


-4.49 




KNCen 


1.5319 


0.867 


13.12'' 


0.926"' 


12.84" 


0.797" 


-4.63 


-5.46 


-7.59 


B 


VWCen 


1.1771 


0.967 


12.80'' 


0.448"' 


12.76" 


0.428" 


-2.93 


-3.85 


-6.08 


B 


XXCen 


1.0395 


0.338 


11.11"' 


0.260"' 


10.90" 


0.266" 


-3.19 


-3.91 


-5.58 


B 


ySDor 


0.9931 


0.529 


7.52'- 


0.040" 


7.50* 


0.052" 


-3.16 


-3.91 


-5.57 




4" Gem 


1.0065 


0.460 


7.78'- 


0.010" 


7.81* 


0.014" 


-3.16 


-3.94 


-5.72 


V 


GHLup 


0.9675 


0.031 


10.05^ 


0.364* 


10.25" 


0.335" 


-2.77 


-3.66 


-5.54 


B 


T Mon 


1.4319 


0.574 


10.82"' 


0.209" 


10.71" 


0.181" 


-4.16 


-5.15 


-7.25 


O 


SMus 


0.9850 


0.266 


9.81'- 


0.147* 


9.57" 


0.212" 


-3.48 


-4.10 


-5.62 





UUMus 


1.0658 


0.865 


12.59'' 


0.413"' 


12.41" 


0.399" 


-3.12 


-3.86 


-5.70 




SNor 


0.9892 


0.343 


9.91"' 


0.189" 


9.87" 


0.179" 


-3.23 


-4.00 


-5.77 


B 


UNor 


1.1019 


0.422 


10.72"' 


0.892"' 


10.46" 


0.862" 


-3.14 


-3.90 


-5.72 




XPup 


1.4143 


0.232 


12.36'' 


0.443* 


11.64" 


0.402" 


-3.57 


-4.38 


-6.34 




APPup 


0.7062 


0.109 


..J 


..J 


...■'' 


../ 


7.37 


6.78 


5.26 


B 


AQPup 


1.4786 


0.436 


12.52"' 


0.512" 


12.41" 


0.518" 


-4.53 


-5.35 


-7.27 


B 


BNPup 


1.1359 


0.397 


12.95'' 


0.438" 


12.93" 


0.416" 


-3.55 


-4.33 


-6.14 




LSPup 


1.1506 


0.012 


13.55'' 


0.478"' 


13.39" 


0.461" 


-3.60 


-4.37 


-6.18 


B 


RSPup 


1.6174 


0.944 


11.56'' 


0.446" 


11.30" 


0.457" 


-4.71 


-5.69 


-7.81 




VZPup 


1.3649 


0.816 


13.08"' 


0.471" 


12.84" 


0.459" 


-3.93 


-4.63 


-6.31 




KQSco 


1.4577 


0.446 


12.36'- 


0.839" 


12.23* 


0.869" 


-4.05 


-5.11 


-7.55 




EUTau 


0.3227 


0.414 


10.27'- 


0.170" 


10.27" 


0.170"' 


-2.26 


-2.74 


-4.05 


Be 


SZTau 


0.4981 


0.744 


8.73" 


0.290" 


8.55" 


0.295" 


-2.38 


-2.93 


-4.33 


B 


TVel 


0.6665 


0.233 


9.80'' 


0.281" 


10.02" 


0.289" 


-2 24 


-2.88 


-4.47 


B 


AX Vel 


0.5650 


0.872 


10.76^ 


0.224* 


...■'' 


0.224" 










RYVel 


1.4496 


0.704 


12.02"' 


0.562" 


11.73" 


0.547" 


-4.23 


-5.05 


-6.96 




RZVel 


1.3096 


0.793 


11.02"' 


0.335" 


10.77" 


0.299" 


-3.78 


-4.61 


-6.56 




SWVel 


1.3700 


0.792 


11.00"' 


0.349" 


11.88" 


0.344" 


-4.02 


-4.83 


-6.75 




SXVel 


0.9800 


0.497 


11A4' 


0.250* 


11.4F 


0.263" 


-3.33 


-3.95 


-5.49 





" Fouque et al. (2007). 

* Benedict et al. (2007). 

■^ Groenewegen et al. (2004). 

"' Storm et al. (2004). 

" Laney & Stobie (1995) and Groenewegen (2004). 

^ Not included in the analysis of the metallicity effect. 

* Groenewegen (2008). 

* Fernie et al. (1995). 



dial velocity correction, which was derived from 20 unblended 
narrow lines well spread over the spectrum, selected among the 
species Fe i, Fe ii and Mg i. The measured signal-to-noise ratios 
vary between 70 and 100 for the FEROS spectra and between 50 
and 70 for the UVES spectra. 

The selected Cepheids span a wide period range, from 3 to 
99 days. For the Galactic Cepheids we have adopted periods, op- 
tical and NIR photometry from Laney & Stobie (1994), Storm 
et al. (2004), Groenewegen et al. (2004), Benedict et al. (2007) 
and Fouque et al. (2007). For the Magellanic Cepheids we have 
adopted periods, optical and NIR photometry from Laney & 
Stobie (1994). The pulsation phases at which our stars were ob- 
served and selected characteristics are listed in Tables 2 and 3 
for the Galactic and the Magellanic Cloud Cepheids, respec- 
tively. Distance and reddening estimates for Galactic Cepheids 
come from two different samples. The "Old Sample" includes 
32 Cepheids (see columns 4 and 5 in Table 2) and among 
them 25 objects have distance moduli provided by Storm et al. 



(2004, see their Table 3) and by Groenewegen et al. (2004, see 
their Table 3). These distances are based on the infrared sur- 
face brightness method and the two different calibrations pro- 
vide, within the errors, the same distances. For the remaining 
seven objects distance estimates are not available in the litera- 
ture. For five of them we determined the distance by combining 
the linear diameter from Laney & Stobie (1995) with the V,K 
unreddened magnitudes from Laney & Stobie (1994), using two 
surface brightness-color calibrations: 

- from Groenewegen (2004): we have combined Eq. 1 and 
Eq. 2 with Table 3 coefficients marked with the filled circle 
(the V vs. V - K relation) 

- from Fouque & Gieren (1997): we have combined Eq. 1 with 
Eq. 27 (the Vvs.V -K relation) 

The distance moduli derived with the two calibrations men- 
tioned above agree very well (within 1 %) and we have adopted 
the distances determined with Groenewegen's calibration (these 
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Table 3. I'ulsation phases (0) and intrinsic parameters of the Magellanic Cepheids. Periods (log P), apparent mean magnitudes 
and reddenings come from Laney & Stobie (1994). The mean /T-band magnitudes were transformed into the 2MASS photometric 
system using the transformation provided by Koen et al. (2007). 



ID 


logP 





Bo 


Vq 


Ko 


E(B - V) 


LMC 


HV877 


1.654 


0.682 


14.06 


12.98 


10.77 


0.12 


HV879 


1.566 


0.256 


14.12 


13.15 


11.03 


0.06 


HV971 


0.968 


0.237 


14.86 


14.24 


12.68 


0.06 


HV997 


1.119 


0.130 


14.94 


14.19 


12.37 


0.10 


HV 1013 


1.382 


0.710 


14.39 


13.46 


1141 


0.11 


HV 1023 


1.425 


0.144 


14.48 


13.51 


11.45 


0.07 


HV 2260 


1.112 


0.144 


15.19 


14.43 


12.67 


0.13 


HV 2294 


1.563 


0.605 


13.19 


12.45 


10.74 


0.07 


HV 2337 


0.837 


0.861 






13.27 


0.07 


HV 2352 


1.134 


0.201 


14.49 


13.84 


12.25 


0.10 


HV 2369 


1.684 


0.136 


13.15 


12.29 


10.38 


0.10 


HV 2405 


0.840 


0.037 






13.43 


0.07 


HV 2580 


1.228 


0.119 


14.33 


13.67 


11.92 


0.09 


HV 2733 


0.941 


0411 


14.85 


14.34 


13.00 


0.11 


HV 2793 


1.283 


0.917 


14.49 


13.58 


11.75 


0.10 


HV 2827 


1.897 


0.880 


13.19 


12.03 


9.80 


0.08 


HV 2836 


1.244 


0.059 


14.85 


14.02 


12.04 


0.18 


HV 2864 


1.041 


0.055 


15.16 


14.42 


12.77 


0.07 


HV 5497 


1.997 


0.321 


12.73 


11.63 


9.43 


0.10 


HV 6093 


0.680 


0.024 


15.74 


15.16 


13.71 


0.06 


HV 12452 


0.941 


0.860 


15.25 


14.60 


12.83 


0.06 


HV 12700 


0.911 


0.342 


15.62 


14.87 


13.12 


-0.01 


SMC 


HV817 


1.277 


0.298 


14.13 


13.59 


12.12 


0.08 


HV823 


1.504 


0.873 


14.46 


13.60 


11.58 


0.05 


HV824 


1.818 


0.315 


13.06 


12.27 


10.33 


0.03 


HV829 


1.931 


0.348 


12.61 


11.81 


9.92 


0.03 


HV834 


1.866 


0.557 


12.95 


12.14 


10.20 


0.02 


HV837 


1.631 


0.822 


13.95 


13.10 


11.11 


0.04 


HV847 


1.433 


0.500 


14.40 


13.66 


11.83 


0.08 


HV865 


1.523 


0.108 


13.55 


12.93 


11.21 


0.06 


HV 1365 


1.094 


0.184 


15.39 


14.79 


13.20 


0.07 


HV 1954 


1.223 


0.847 


14.13 


13.62 


12.12 


0.07 


HV 2064 


1.527 


0.279 


14.28 


13.50 


11.61 


0.07 


HV 2195 


1.621 


0.135 


13.85 


13.07 


11.09 


-0.02 


HV 2209 


1.355 


0.822 


13.99 


13.42 


11.84 


0.04 


HV 11211 


1.330 


0.516 


14.36 


13.64 


11.83 


0.06 



are the values listed in Table 2). The "New Sample" includes 
32 Cepheids (see columns 6 and 7 in Table 2) and among them 
24 objects have distance moduli based on the infrared surface 
brightness method provided by Fouque et al. (2007, see their 
Table 7). The trigonometric parallaxes for 1 Car, /3 Dor, and 
^ Gem have been provided by Benedict et al. (2007). For EU 
Tau, we adopted the distance by Groenewegen et al. (2004), 
while for KQ Sco and SX Vel we adopted distances calculated 
by one of us (MG), following the general method outlined in 
Groenewegen (2007), but using the SB-relation and projection- 
factor from Fouque et al. (2007) for consistency. The Cepheid 
AP Pup was not included in the analysis of the metallicity ef- 
fect because an accurate estimate of its linear diameter is not 
available. For this object in Table 2 we only listed the apparent 
magnitudes. The same outcome applies to AX Vel, since Laney 
& Stobie (1995) mentioned that the quality of the radius solu- 
tion for this object was quite poor. Accurate reddening estimates 
for Galactic Cepheids have been recently provided by Laney & 
Stobie (2007), however, we typically adopted the reddening es- 
timates used to determine individual Cepheid distances. 

The current Cepheid sample includes objects that are clas- 
sified in the literature as fundamental pulsators, except EU Tau 



and SZ Tau, which are classified as first overtone pulsators, and 
AX Vel that is classified as one of the few double-mode pul- 
sators in the Galaxy (Fernie et al. 1995; Sziladi et al. 2007). 
Indeed, when plotting our Cepheid sample in the log P vs My 
plane, we confirm that EU Tau and SZ Tau are the only stars 
lying on the first overtone PL relation. Therefore, their ob- 
served periods have been "fundamentalised" using the relation 
Pq = Pi/(0.716 - 0.027 logPi) (Feast & Catchpole 1997). 



Approximately, 50% of Galactic Cepheids in our sample are 
spectroscopic binaries (Szabados 2003) and one star (^ Gem) is a 
visual binary (see last column of Table 2). According to current 
empirical evidence the companions are typically B and A-type 
main sequence stars, which are much less luminous (at least 3 
mag) than our main Cepheid targets. Only in the case of the two 
most luminous B dwarf, companions respectively of KN Cen 
and S Mus, we have detected a small contribution to the con- 
tinuum level of the Cepheid spectra. Its effect on the final iron 
abundances will be discussed in Section 4. 
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Fig. 2. Examples of visual inspection on the selected iron lines. 
The iron lines are plotted in black while the other elemental lines 
are plotted in gray. In the top panel it is shown an unblended line, 
while in the bottom panel there is an example of a blended line. 



3. Methodology 

Below we present the methodology used to assemble the linelist 
and to determine the equivalent widths and the atmospheric pa- 
rameters of our programme stars. All of these are key ingredients 
in the calculation of elemental abundances. 



3.1. Linelist 

A crucial step of any spectral analysis in order to derive elemen- 
tal abundances is a careful assembling of a line list. We have 
assembled our Fei - Fen line list from Clementini et al.(1995), 
FC97, Kiss & Vinko (2000) and Andrievsky et al.(2002a) plus 
a selection of lines from VALD (Vienna Atomic Line Database; 
Kupka et al. 1999; Ryabchicova et al. 1999).The VALD lines 
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Fig. 3. The distributions of the equivalent widths (EW) of our 
lines (empty histogram) and Fry & Carney's lines (1997, dark 
grey histogram), in the case of the same star (VCen). 



have been selected for effective temperatures typical of Cepheid 
stars (4500-6500 K). We have, then, visually inspected each line 
on the observed spectra, in order to check their profile and to 
discard blended lines. In order to do so, we have searched the 
VALD database (with the command extract stellar) for all the 
existing lines between 4800 and 7900 A on stellar spectra char- 
acterized by stellar parameters typical of Cepheids (Teff'=4500 
K, 5500 K and 6500 K, log g=2 and vt=3km /s). We have, then, 
overplotted all the lines found in VALD, that fall within ± 3 A 
from each of our iron lines, and checked for their possible con- 
tribution to the equivalent width of our iron line (an example 
is shown in Fig. 2). This test was carried out on 3 different ob- 
served spectra, characterized by effective temperatures (as found 
in the literature from previous analyses) close to 4500 K, 5500 K 
and 6500 K, respectively. Any contribution larger than 5% of the 
line strength made us discard the iron line under scrutiny. There 
remains, of course, the possibility that some weak lines could be 
missing in the VALD compilation, but these are expected to be 
only weak lines, whose contribution would, then, be negligible. 

Our final list includes 275 Fe i lines and 37 Fe ii lines, span- 
ning the spectral range covered by our FEROS spectra. Our 
UVES spectra cover a narrower spectral range for which we can 
use 217 Fei lines and 30 Feii lines. For all the lines, we have 
adopted the physical properties (oscillator strengths, excitation 
potentials) listed in VALD. Figure [3] shows the comparison be- 
tween the distribution of the equivalent widths of our lines with 
those from FC97 for the Galactic Cepheid V Cen. As it can be 
clearly seen, our list is significantly larger and well samples the 
best range of equivalent widths, around 20 mA, as suggested by 
Cayrel (1988) to obtain reliable elemental abundances (iron in 
our case). The line list is presented in Table 4, where the four 
columns list respectively the wavelength, the ion identification, 
the excitation potential and the gj- values of each line. 
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Table 4. Iron line list, from left to right the columns display 
wavelength, ion identification, excitation potential (EP) and log 
Rvalues. The table is available in its entirety via the link to the 
machine-readable version. A portion is shown here for guidance 
regarding its form and content. 



A 


Ion 


EP 


logg/ 


4892.87 


Fei 


4.22 


-1.34 


4917.23 


Fei 


4.19 


-1.29 


4924.77 


Fei 


2.28 


-2.24 


4932.08 


Fei 


4.65 


-1.49 


4950.10 


Fei 


3.42 


-1.67 


4973.11 


Fei 


3.96 


-0.95 


4994.13 


Fei 


0.91 


-3.08 


5029.62 


Fei 


3.41 


-2.05 


5044.21 


Fei 


2.85 


-2.04 


5049.82 


Fei 


2.28 


-1.36 


5054.64 


Fei 


3.64 


-1.92 


5056.84 


Fei 


4.26 


-1.96 


5067.15 


Fei 


4.22 


-0.97 


5072.67 


Fei 


4.22 


-0.84 


5083.34 


Fei 


0.96 


-2.96 



quality of the data (Equation 7, Cayrel 1988). We can then 
safely use all our EW measurements, independently of the 
method we used to derive them. 

For the determination of the metallicities, we have selected 
only lines with equivalent widths between 5 and 150 mA. The 
lower limit was chosen to be a fair compromise between the 
spectral characteristics and the need for weak lines for an op- 
timal abundance determination. The upper limit was selected in 
order to avoid the saturated portion of the curve of growth. We 
note that a slightly higher upper limit (200 mA) was chosen for 
T Mon, SZ Tau and KQ Sco because these stars have very strong 
lines. This was done in order to keep the number of selected lines 
similar to the one used for the rest of the sample, after checking 
that this higher upper limit does not have any effect on the final 
metallicities derived for these three stars. 

Considering the mean difference on the EW obtained with 
fitline and splot and the uncertainty in the continuum placement 
from two measurements of the EW (carried out independently by 
two of us) we assume + 3 mA as error on the equivalent widths 
for the lines below 100 mA and + 5 mA for the stronger features. 



3.2. Equivalent widths 

The second crucial step of our spectral analysis is the measure- 
ment of the equivalent widths (EW) of the iron lines we have 
assembled in our final list. At first, several independent and man- 
ual (i.e. with the IRAF splot task) measurements of the EW 
of the whole set of Fe i and Fe ii lines in a selected number of 
Cepheids were performed in order to test the reproducibility of 
our measures. However, because of the large number of lines 
and spectra, our final EW measurements were derived by us- 
ing a semi-interactive routine (fitline) developed by one of us 
(PF, fitline). This code is based on genetic algorithms, which 
mimic how DNA can be affected to make the evolution of species 
(Charbonneau 1995). It uses a Gaussian fit, which is defined by 
four parameters: central wavelength, width, depth and contin- 
uum value of the line. A top-level view of the algorithm is as 
follows: 

1 . Compute an initial set of Gaussians, picking random values 
of the four parameter (scaled to vary between and 1) and 
calculate the;i'^ with the observed line for each Gaussian. 

2. Compute a new set of Gaussian from the 20 best fit of the 
previous "generation" introducing random modification in 
the values of the parameters ("mutation'). 

3. Evaluate the "fitness" of the new set (x^ calculation for each 
Gaussian) and replace the old set with the new one. 

4. Iterate the process (100 to 200 "generations") to get the best 
fit (lowest ;\f^) for each observed line. 

All the UVES spectra have been smoothed to improve the 
quality of the Gaussian fit (smooth_step =11 pix using splot 
task in IRAF). All the useful information from the spectra are 
preserved in the process. The selected iron lines, even on the 
smoothed spectra, do not show any contamination from other 
line. However, for about 15% of the lines, the Gaussian profile 
performed by fitline, could not satisfactorily reproduce the 
observed profile. The equivalent widths of these lines, usually 
very broad or asymmetric, were measured manually with the 
splot task. The mean difference, as computed for those lines 
for which both methods could be applied, is around 1.5 mA, 
comparable to the average error on the EW inferred from the 



3.3. Stellar parameters 

The stellar parameters were derived spectroscopically. We have 
determined the stellar effective temperature Teff by using the line 
depth ratios method described in Kovtyukh & Gorlova(2000). It 
is based on weak neutral metal lines (in our case, less than 100 
mA in equivalent width in order to exclude line broadening ef- 
fects) with low excitation potentials, selected to obtain as close 
as possible pairs in wavelength space. Since the calibration of 
this method has been done using the FC97 scale of tempera- 
tures, our temperature scale is linked to theirs. The line depth ra- 
tios have the advantage of being independent of interstellar red- 
dening and metallicity effects, uncertainties that instead plague 
other methods like the integrated flux method or colors calibra- 
tions (Gray 1994). The main uncertainties of the calibration of 
the line depth ratios, instead, lie in the accuracy of its zero point 
and slope, which have been thoroughly tested either with differ- 
ent colour-temperature relations or diameters measurements by 
Kovtyukh & Gorlova (2000). As a sanity check on our tempera- 
ture determinations, we have analysed five individual spectra of 
Galactic Cepheid AP Puppis taken at different phases along the 
pulsation period (0=0.11 - observed twice -, 0.31, 0.51, 0.91). 
Four out of these five spectra have been collected in an indepen- 
dent observational campaign and the results have already been 
discussed by Lemasle et al. (2007). Figure |4] shows the effec- 
tive temperatures derived for the five phases of AP Pup (filled 
squares) in comparison with the temperature variation estimated 
by Pel (1978) on the basis of accurate Wakaven photometry. The 
agreement is indeed very good. 

The total number of line depth ratios adopted to estimate the 
temperature ranges from 26 to 32 and from 20 to 28 for Galactic 
and Magellanic Cepheids, respectively. From this method, we 
have obtained effective temperatures with an intrinsic accuracy 
of about 40 K for the Galactic Cepheids and 50 K for the 
Magellanic ones (errors on the mean). Table 7 and Table 9 
list our final effective temperatures for Galactic and Magellanic 
Cepheids, respectively. 

Microturbulent velocity (v,) and gravity (log^) were con- 
strained by minimizing the log([Fe/H]) vs EW slope (using the 
Fe I abundance) and by imposing the ionisation balance, respec- 
tively. These two procedures are connected and require an it- 
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Fig. 4. Comparison between the behaviour of the effective tem- 
perature of the Galactic Cepheid AP Puppis found by Pel (1978, 
solid line) and our results using the line depth ratios method for 
5 phases (filled squares) with their associate errorbars. 



erative process (on average, between 5 and 7 iterations, de- 
pending on the star). We first achieved the minimization of the 
log([Fe/H]) vs EW slope and, subsequently, the ionisation bal- 
ance. As first guess for the microturbulent velocity and the grav- 
ity, we adopted values typical of Cepheids (v,= 3 km/s, \ogg - 
2) as inferred from previous studies (FC97 and Andrievsky et al. 
2002a, 2002b). 

We first assumed the ionisation balance to be fulfilled when 
the difference between [Fe i/H] and [Fe ii/H] is less than the stan- 
dard error on [Fe ii/H] (typically, 0.08-0. 1 dex). However, we no- 
ticed that for most stars this condition was usually satisfied by 
more than one value of log g, suggesting that our conditions for 
fulfillment might be too conservative and not very informative. 
We then checked which log g value satisfies the ionisation bal- 
ance within the standard error on [Fei/H] (typically ,0.02 dex). 
For 55 stars out of 68 we were able to reach the Fe i-Fe ii bal- 
ance within few hundredths of a dex. The corresponding log g 
values are the final gravity values quoted in Table 7 and Table 
9. For the remaining 13 objects, we assumed as our final \ogg 
the value giving the "best" ionisation balance, i.e. the one with 
the smallest difference between [Fei/H] and [Feii/H]. Also, it 
is worth noticing that for two stars (HV 2827and HV 11211) 
we had to increase the temperatures as determined from the line 
depth ratios by 50 K (which is still within the estimated error 
on the determination of Teff), in order to fulfill our requirements 
for a satisfactory ionisation balance. Moreover, we note that the 
log g values derived for the five different phases of AP Puppis 
follow well the trend found by Pel (1978), as it can be seen in 
Fig. |5] At 0=0.11, we have obtained the same value from the 
analysis of the two spectra. 

In order to determine the errors on the microtutbulent veloc- 
ity and the gravity, we have run several iterations for each star, 
slightly modifying the values of these two quantities that fulfill 
the requirements mentioned above. We have estimated errors in 
V, to be 0. 1 km/s and in log g to be 0. 10 dex. 
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Fig. 5. Comparison between the behaviour of the gravity of the 
Galactic Cepheid AP Puppis found by Pel (1978, the two solid 
lines indicate the upper and lower limits) and our results for 5 
phases (filled squares) with their associate errorbars. At 0=0.1 1, 
we have obtained the same value from the analysis of the two 
spectra. 



3.4. Model atmospheres 

We have derived the iron abundances of our stars by using the 
Kurucz WIDTH9 code (Kurucz 1993) and LTE model atmo- 
spheres with the new opacity distribution functions (ODFs) com- 
puted by Castelli & Kurucz (2003). These models neglect en- 
velope overshooting. We have used a grid of solar metallicity 
models for the Galactic and LMC Cepheids and a grid of mod- 
els computed assuming [Fe/H]=-1.0 and an a-element enhance- 
ment of +0.4 dex for the SMC Cepheids (Gratton, Sneden & 
Carretta, 2004). The grids of models have been interpolated in 
temperature in order to match the effective temperature derived 
for each star with the line depth ratio method and in 0.10 dex 
gravity steps. 

Our choices of model atmospheres (in terms of ODFs, over- 
shooting, metallicity, and or-enhancement) have been thoroughly 
tested. A comparison between model atmospheres computed 
with the new ODFs (Castelli & Kurucz 2003, hereafter identi- 
fied as the 2003 models) and the old ODFs (Castelli et al 1997, 
hereafter identified as the 1997 models) shows small differences 
in the derived iron abundances: ~ 0.01-0.03 dex. We tested the 
effect of different treatments of overshooting as implemented in 
different versions of Kurucz models by running the WIDTH9 
code with two sets of the 1997 models computed with and with- 
out the "approximate" overshooting. Differences in the final iron 
abundances amount to 0.05 dex for Galactic and SMC Cepheids 
and around 0.06 dex for LMC Cepheids, with the models with- 
out overshooting giving the lower abundances. 

Since the LMC has a mean metallicity around -0.3 dex,for 
this galaxy we also tested model atmospheres computed for 
[Fe/H]=-0.5 dex, finding differences of the order of ~ 0.01-0.02 
dex in the derived iron abundances compared to the solar metal- 
licity models. Similar differences are found when the [Fe/H]=- 
0.5 dex models are used also for the Small Magellanic Cloud (its 
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Table 5. Effects on measured Fe i and Fe ii abundances caused 
by changes in atmospheric parameters. 




Fig. 6. Histograms of the [Fe/H] ratios derived for all the stars 
of our sample in the Galaxy, the LMC and the SMC. The solid 
lines indicate the mean values found using F and K supergiants 
by Hill et al. (1995) and Hill (1997) (see Table 8 and Section 
4.3). We found a good agreement with the results obtained for 
non-pulsating stars of the same age of Cepheids. 



mean metallicity is around -0.7 /-O .8 dex). No differences in the 
derived iron content have been found between models computed 
with and without the +0.4 dex a-enhancement. 

In all our computations, we have adopted a solar iron abun- 
dance of log[n(Fe)] = 7.51 on a scale where log[n(H)] - 12 
(Grevesse & Sauval, 1998) and we have assumed our final stel- 
lar metallicity to be the Fe i value, which has been derived by a 
far larger number of lines with respect to Fe ii. 



A[Fei/H] A[Feii/H] 



Areff=+100K 

Alogg =+0.1 dex 

Av,=+0.1 Km/s 



+0.07 dex 
+0.00 dex 
-0.03 dex 



+0.00 dex 
+0.04 dex 
-0.02 dex 



ered by our spectra, these bright companions give a contribution 
only to the continuum level, because all their absorption lines 
fall in the ultra-violet region (these are very hot stars, with ef- 
fective temperature ~ 20,000 K). In order to test the effect of 
this contribution, we have subtracted from the observed spec- 
trum of the two binaries the estimated optical spectrum of the 
B dwarf, which only consists of a continuum without any line. 
On the resulting spectra, which we assume to be the true spectra 
of our Cepheids, we have then remeasured the EWs of a sub- 
set of iron lines and derived the metallicity. In the case of KN 
Cen the differences in the EWs measured on the observed spec- 
trum and the true Cepheid spectrum are negligible (i.e. within 
the errors), thus the iron content we have derived for this star is 
robust. Regarding S Mus, instead, the EWs measured on the true 
spectrum are 12-15% larger than the ones measured on the spec- 
trum we have observed. Thus, the true iron abundance of this 
Cepheid star should be 0. 1 dex higher than the one we have de- 
rived. This happens because the intensity of the lines (due only 
to the Cepheid contribution) when compared to the continuum 
of the combined spectra is, in percentage, less than the intensity 
of the lines compared to the continuum of the Cepheid alone. In 
other words the contribution of the companion makes the lines 
of the Cepheid weaker. Please note that the iron content reported 
for S Mus in Table 7 has already been corrected for the above 
mentioned effect. 



4. Iron abundances 

Our final iron abundances, together with the adopted stellar pa- 
rameters, are listed in Table 7 and Table 9 for the Galactic and 
Magellanic Cepheids, respectively. 

For our Galactic sample, we find that the mean value of the 
iron content is solar (0 = 0. 10, see Fig. 6), with a range of values 
between -0.18 dex and +0.25 dex. These two extremes are repre- 
sented respectively by 2 and 1 stars (out of the 32 in total). Our 
Galactic Cepheids span a narrow range of Galactocentric dis- 
tances (the mean Galactocentric distance of our sample is 7.83 
+ 0.88 kpc, see Fig. 7), thus preventing us from giving any indi- 
cations about the metallicity gradient in the Galactic disc. 

For the LMC sample, we find a mean metallicity value of ~ 
-0.33 dex ((f) = 0.13, see Fig. 6), with a range of values between - 
0.62 dex and -0. 10 dex. Here, the more metal-rich extreme is just 
an isolated case, while the metal-poor end of the distribution is 
represented by 3 stars. 

For the SMC sample, we find that the mean value is about ~ 
-0.75 dex (0 - 0.08, see Fig. 6), with a range of values between 
-0.87 and -0.63. 

As we have already mentioned in Section 2, there are two 
binary stars (KN Cen and S Mus) among our Galactic sample 
with a bright B dwarf as a companion. In the spectral range cov- 



Rg [Kpc] 

Fig. 7. Histograms of the Galactocentric distances derived for all 
the Galactic Cepheids with individual heliocentric distances. 
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Table 6. Stellar parameters and iron abundances along the pul- 
sation cycle of the Galactic Cepheid AP Puppis. Phase 0. 11 was 
observed twice 50 days apart, leading to very consistent results. 



Phase 


Teff 


V, 


log(g) 


[Feim 


0.11 


6070 


3.05 


2.1 


-0.07 


0.11 


6130 


3.00 


2.1 


-0.03 


0.31 


5750 


2.80 


1.9 


-0.03 


0.51 


5550 


3.30 


2.0 


-0.06 


0.91 


6160 


4.20 


2.2 


-0.13 



4.1. Uncertainties 

The internal uncertainties in the resulting abundances are due 
to errors in the atomic data (^/-values) and EW measurements. 
We have estimated, on average, an internal error in the [Fe/H] 
determination of 0. 10 dex. 

It is also important to understand the effects of potential 
systematic errors in the stellar parameters on the final derived 
abundances (see Table 5). In order to do so we have determined 
curves of growth for different effective temperatures, microtur- 
bulent velocities and gravities for both Fei and Fen. As ex- 
pected, we find that Fe i abundances marginally depend on the 
gravity, whereas Fe ii does not depend on the effective tempera- 
ture. An increase in temperature of 100 K, at fixed v, and logg, 
results in an increase in [Fe i/H] of about 0.07 dex. An increase 
in V, of 0. 1 km/s gives a decrease in [Fe i/H] of about 0.03 dex, 
for constant Tf^g and logg, and we obtain a decrease of 0.02 dex 
for [Feii/H]. An increase in logg of about 0.1 dex produces an 
increase in [Fe ii/H] of about 0.04 dex (with fixed v, and Tetf). 

An additional potential source of uncertainty and concern 
comes from the fact that our metallicities have been derived from 
single epoch observations. However, we note that FC97 did not 
find any significant difference in their derived [Fe/H] as a func- 
tion of phase (the test was performed on four of their longest 
period cepheids) and Luck & Andrievsky (2004) and Kovtyukh 
et al. (2005) show that the elemental abundances for Cepheids 
with a period between 6 and 68 days are consistent for all pul- 
sation phases. Moreover, our exercise on AP Puppis (see section 
3.3 and Table 6) further confirms these conclusions, for a star 
even with a shorter period (5 days). As it can be seen in Table 
6, not only do the derived stellar parameters nicely follow the 
expectations from photometry, as shown in Fig. |4] and Fig.|5] but 
the metallicities measured at all phases agree within the errors 
(-0.1 dex). 

4.2. Comparison witli previous studies on Ceplieids 

Before we address the key issue concerning the metallicity de- 
pendence of the Cepheid PL relation (cf. Sect. 5), let us first 
compare our results to previous works. For this purpose, we 
have selected the chemical analyses of FC97, Andrievsky et al. 
(2002a, 2002b, 2002c) and Luck et al. (2003) for the Galactic 
Cepheids and LL92 and L98 for the Magellanic Clouds. When 
necessary, we have rescaled the literature [Fe/H] values to the 
solar iron abundance we have adopted {log[n{Fe)] - 7.51). One 
should keep in mind that the same stars have likely been ob- 
served at different phases and analysed with different tools and 
model atmospheres, which may lead to the determination of dif- 
ferent combinations of stellar parameters. In the case of obser- 
vations at different phases, as we already mentioned in the pre- 
vious section, our results on the five spectra of AP Pup and the 
conclusions of Luck & Andrievsky (2004) and Kovtyukh et al. 



(2005) show that the elemental abundances do not depend on 
the phase. When multi-phase observations of the same star were 
available in the literature (this is the case for some of the Galactic 
Cepheids), we have compared the [Fe/H] values obtained for the 
same phase of our observations. Otherwise, we chose the iron 
content corresponding to the closest value of effective temper- 
ature to ours. Regarding the determination of different sets of 
stellar parameters, we do not have the capabilities (in terms of 
analytical codes and different sets of model atmospheres) in or- 
der to properly take this into account. However, we did signal 
such differences whenever we noted them. 

In total, for the Galactic sample, we have 6 stars in common 
with FC97 and 18 with the entire sample of Andrievsky 's group, 
while for the MC sample, we have 3 stars in common with LL92 
and 7 with L98. 



4.2.1. Galactic Cepheids 

The mean difference between our results and those of FC97 and 
Andrievsky analyses is comparable to the difference between 
FC97 and Andrievsky's values (0.08 + 0.02). In more detail, the 
mean difference between our iron abundances and FC97 is 0.09 
+ 0.02, which is satisfactory. For 4 stars (V Cen, S Nor, TMon, 
( Gem) the agreement is at Icr level, for the remaining 2 stars 
(we have 6 in common, in total) the agreement is well within 
the quoted uncertainties. We note that our derived abundances 
for V Cen and S Nor are more in agreement with the metallicity 
derived by Andrievsky et al. than with FC97. 

When comparing our results to Andrievsky's, we obtain a 
mean difference of 0.07 + 0.05. For 14 stars (out of 18), the 
iron abundances agree quite well within the associated errors. 
Of the remaining 4 stars, we note that SW Vel, /3 Dor, and T 
Mon agree within the standard deviation (cr=0.20, 0. 19, 0.22 dex 
respectively) instead of the standard error, that is the condition 
suggested by Kovtyukh et al. (2005a) for the independence of 
the elemental abundance on the phase. This is not the case for ( 
Gem, for which we cannot explain the difference (0.20 dex) but 
for which we find instead an agreement with FC97. 

4.2.2. Magellanic Ceplieids 

In order to properly compare our results for the Magellanic 
Cepheids with the values obtained by LL92 and L98, we first 
had to rescale the latter values for the difference in the adopted 
solar iron abundance between us (7.51) and them (7.67 and 7.61, 
respectively). These revised values are listed in the last column 
of Table 9. 

In general, the mean metallicities that L98 found with their 
complete sample (-0.30 dex and -0.74 dex for the LMC and 
SMC, respectively) are in very good agreement with our results 
(-0.33 dex and -0.75 dex). They also found a similar spread in 
iron for the MCs. With LL92, instead, there is a good agree- 
ment in the case of LMC but a difference for the metallicity of 
the SMC (0.15 dex greater than our mean value). Our derived 
abundances are always smaller than the values derived by LL92. 
However, the number of objects in common is too small to con- 
strain the effect on a quantitative basis. 

When we move to a star-by-star comparison, larger differ- 
ences emerge. The comparison with L98 abundances discloses 
a very good agreement for one object (HV 837) and a plausi- 
ble agreement, i.e. within the standard deviation error, for other 
3 stars (HV 5497, HV 824 and HV 834). Regrettably, this is 
not the case of HV 879, HV 2827 and HV 829 for which we 
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Table 7. Stellar parameters and Fe i and Fe ii abundances of our Galactic sample. The Fe i values have been adopted as final [Fe/H]. 
We compare our results, where it is possible, with previous investigations: Fry & Carney (1997); Andrievsky et al. (2002a,b,c), Luck 
et al. (2003). 



ID 


Teff 


Vf 


logs 


[Fei/H] 


[Feii/H] 


[Fe/H]fc 


[Fe/H]^ 


ICar 


4750 


3.60 


0.4 


0.00 


0.00 






UCar 


5980 


3.00 


1.2 


+0.17 


+0.13 






VCar 


5560 


3.05 


1.8 


-0.07 


-0.02 






WZCar 


4520 


1.80 


2.1 


+0.08 


+0.19 






VCen 


6130 


2.80 


1.9 


+0.04 


-0.03 


-0.12 


+0.04" 


KNCen 


5990 


3.80 


1.6 


+0.17 


+0.09 






VWCen 


5240 


4.20 


1.2 


-0.02 


+0.07 






XXCen 


5260 


2.95 


1.3 


+0.04 


+0.04 






jSDor 


5180 


3.00 


1.3 


-0.14 


-0.11 




-O.or 


^Gem 


5180 


3.70 


1.4 


-0.19 


-0.14 


-0.04 


+0.06° 


GHLup 


5480 


3.60 


1.5 


+0.05 


+0.01 






TMon 


4760 


3.80 


0.6 


-0.04 


-0.03 


+0.09 


+0.21* 


SMus 


5780 


2.75 


1.8 


+0.13 


+0.19 






UUMus 


5900 


3.05 


1.7 


+0.11 


+0.05 






SNor 


5280 


2.80 


1.1 


+0.02 


+0.04 


-0.05 


+0.06" 


UNor 


5230 


2.60 


1.1 


+0.07 


+0.13 






XPup 


5850 


3.30 


1.4 


+0.04 


-0.05 




0.00" 


APPup 


6070 


3.05 


2.1 


-0.07 


-0.07 






AQPup 


5170 


3.10 


0.8 


-0.07 


-0.09 




-0.14* 


BNPup 


5050 


2.95 


0.6 


-0.10 


-0.07 




+0.0L 


LSPup 


6550 


3.50 


2.2 


-0.16 


-0.10 






RSPup 


4960 


3.50 


0.7 


+0.09 


+0.10 




+0.16* 


VZPup 


5230 


3.25 


1.1 


-0.17 


-0.15 




-0.16" 


KQSco 


4840 


3.60 


0.7 


+0.21 


+0.27 




+0.16'' 


EUTau 


6060 


2.30 


2.1 


+0.04 


+0.02 




-0.03" 


SZTau 


5880 


2.80 


1.7 


+0.07 


+0.04 


-0.02 


+0.06" 


TVel 


5830 


2.55 


1.8 


+0.10 


+0.03 




-0.04* 


AXVel 


5860 


3.10 


1.8 


+0.02 


-0.06 






RYVel 


5250 


4.10 


1.2 


-0.07 


-0.06 




-0.03* 


RZVel 


5140 


4.40 


1.6 


-0.19 


-0.14 




-0.11* 


SWVel 


6590 


3.75 


1.7 


-0.24 


-0.17 


-0.10 


-0.05* 


SXVel 


5380 


2.55 


1.2 


+0.02 


-0.02 




-0.04* 


" Luck et al. (2003). 












* Andrievsky et al. 


(2002a). 












^ Andrievsky et al. 


(2002c). 













Andrievsky et al. (2002b). 



note rather large discrepancies (~ 0.2 - -0.3 dex) that remain un- 
explained. However, we note that L98 used different analytical 
codes, different oscillator strengths, and different values of the 
stellar gravity (they adopted the "physical" gravity calculated 
from the stellar mass, the luminosity and the temperature). A 
combination of all these three factors could well account for the 
observed differences. 

As already mentioned, larger differences (0.4 - 0.6 dex) are 
found with respect to LL92, for the 3 stars we have in common. 
Following the referee suggestion we performed a more quantita- 
tive comparison for the three Cepheids in common with LL92, 
namely HV 865, HV 2195 and HV 2294. The calibrated spec- 
tra were kindly provided in electronic form by R.E. Luck. By 
adopting the same atmospheric parameters used by LL92 we 
find the following iron abundances: [Fei/H]= -0.75 (HV 865), 
-0.60 (HV 2195) and -0.15 (HV 2294). The comparison with 
the iron abundances provided by LL92 (see the last column in 
Table 9) indicates that the new measurements are ~ 0.3 dex 
systematically more metal-poor. They also agree, within the er- 
rors, reasonably well with current measurements (see column 5 
in Table 9). Thus supporting the arguments quoted above to ac- 
count for the difference between the two different iron measure- 
ments. Moreover, one should keep in mind that the quality of 



the LL92 data set is significantly lower than ours (R x 18,000 
vs 40,000) and that their oscillator strengths may differ from 
our selection: they were taken from the critical compilations of 
Martin, Fuhr & Wiese (1988) and Fuhr, Martin & Wiese (1988), 
which are included in VALD but not among the references used 
for our gf-values. Moreover, they have used different analytical 
tools: MARCS model atmospheres (Gustafsson et al. 1975) and 
a modified version of the LINES line-analysis code and MOOG 
(Sneden 1973) synthesis code. They also state that their method 
overestimates the equivalent widths of the weak lines. 



4.3. A comparison with different stellar populations in the 
Magellanic Clouds 

The mean metallicities of our Magellanic sample are in good 
agreement with previous results obtained for F, K supergiants 
and B stars in the Magellanic Clouds (see Table 8). 

Russell & Bessell (1989) found a mean [Fe/H]= -0.30 + 
0.2dex in the LMC and [Fe/H]= -0.65 ± 0.2 dex in the SMC, 
analysing high-resolution spectra of 16 F supergiants (8 for each 
galaxy). In 1995 Hill et al. obtained, from 9 F supergiants from 
the field of the LMC, a mean iron abundance of -0.27 + .15 dex 
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Table 8. Comparison of the mean metallicities of the Magellanic 
Clouds with previous studies. The number of studied stars are 
also listed. RB89: Russell & Bessell (1989), RD92: Russell 
& Dopita(1992), R93: Rolleston et al. (1993), H95: Hill et 
al. (1995), H97: Hill (1997), KOO: Korn et al. (2000), AOl 
Andrievsky et al. (2001), R02: Rolleston et al. (2002), G06 
Grocholski et al. (2006), T07: Trundle et al. (2007), P08 
Pompeia et al. (2008). 



Reference 


[Fe/H]i„c 


[Fe/HlsMC 


Notes 


This worlc 


-0.33+0.13 


-0.75 + 0.08 


22+14 Cepheids 


RB89 


-0.30±0.20 


-0.65 + 0.20 


8+8 F supergiants 


R93 




-0.80 + 0.20" 


4 B stars 


H95 


-0.27+0.15 




9 F supergiants 


H97 




-0.69 + 0.10 


6+3 K supergiants 


KOO 


-0.40+0.20 


-0.70 + 0.20 


6 B stars 


AOl 


-0.40+0.15 




9 F supergiants 


R02 


-0.31+0.04 




5 B stars 


G06 


-0.3 - -2.0* 




200 giants 


TO? 


-0.29+0.08 


-0.62+0.14 


13+5 cluster giants 


P08 


-0.3- -1.8 




59 red giants 











" Mean metallicity based on oxygen. 
* Metallicity based on calcium triplet. 



and Hill (1997) found a mean [Fe/H]= -0.69 + 0.1 dex, analysing 
6 K supergiants in the SMC. Andrievsky et al. (2001) re-analysed 
the sample of F LMC supergiants studied by Hill et al. (1995) 
and obtained a slightly lower mean value: [Fe/H]= -0.40 + 0.15 
dex. 

Regarding the B type stars, Rolleston et al. (1993) and 
Rolleston et al. (2002) found, from the analysis of 4 stars in 
the SMC and 5 stars in the LMC, respectively, mean values of 
metallicity of -0.8 + 0.20 dex and -0.31 + 0.04 dex. Korn et 
al.(2000) obtained for the LMC (from 6 B stars) a mean [Fe/H] = 
-0.40 + 0.2 dex and [Fe/H]=-0.70±0.2 dex for the SMC (from 
3 B stars). Grocholski et al. (2006) using homogeneous cal- 
cium triplet measurements for 200 stars belonging to 28 differ- 
ent LMC clusters covering a broad range of cluster ages, found 
by transforming the Ca abundance in iron abundance that the 
metallicity ranges from [Fe/H] — 0.3 to [Fe/H]~-2.0. More re- 
cently. Trundle et al. (2007) collected high resolution spectra 
with FLAMES @VLT for a good sample of Magellanic giants 
and found a mean metallicity of [Fe/H=-0.62 + 0.14, 5 giants 
in the SMC cluster NGC 330, and of [Fe/H] =-0.29 ± 0.08, 13 
giants in the LMC cluster NGC 2004. By using similar data for 
59 red giant stars in the LMC inner disk, Pompeia et al. (2008) 
found that their metallicity ranges from [Fe/H]~0.0 to [Fe/H]~- 
1.80. 

In conclusion, the mean iron content of Cepheids, for the 
Magellanic Clouds, agree very well with the results obtained for 
non variable stars of similar age and stars that are Cepheid's pro- 
genitors. Cepheids do not show any difference with these two 
other populations. 

5. The effect of [Fe/H] on the PL relation 

To assess the effect of the iron content on the Cepheid PL rela- 
tion we select, among our sample, only the stars with periods be- 
tween 3 and 70 days (61 stars out of 68), populating, in this way, 
the linear part of the PL relation, i.e. the one useful for distance 
determinations (e.g. Bono et al. 1999; Marconi et al. 2005). 

By using the values given in the current literature (e.g. 
Benedict et al. 2002; Walker 2003; Borissova et al. 2004; 



SoUima et al. 2008; Catelan & Cortes 2008; Groenewegen et 
al. 2008), we adopted for the barycentre of the LMC a true dis- 
tance modulus (p-LMc) of 18.50 mag (with an error of ± 0.10). 
This is also consistent with the standard PL relations (Freedman 
et al. 2001; Persson et al. 2004) used as comparison. The SMC 
is considered 0.44 + 0.10 mag more distant than the LMC (e.g. 
Cioni et al. 2000; Bono et al. 2008). This value of the relative 
distance between the two galaxies has confirmed the results of 
previous studies (Westerlund 1997, and references therein). For 
all the Galactic Cepheids we have individual distances, as men- 
tioned in Section 2 above. 

It is necessary to divide our Cepheid sample into bins of 
metallicity to investigate its effect on the PL relation. The num- 
ber of bins needs to be chosen taking into account two compet- 
ing effects. On the one hand, dividing the stars in more bins in 
principle allows to disentangle finer details. On the other, how- 
ever, each bin needs to contain enough stars so that the instability 
strip is well populated and, hence, spurious sampling effects are 
negligible. We choose two bins with about 30 stars each as best 
compromise between a detailed investigation and statistical sig- 
nificance. This choice is justified by the simulations described in 
the next section. 

Our results for the V and K bands are presented in Figures 
8 and 9, respectively. In the top panels are shown the PL rela- 
tions in each metallicity bin calculated with a linear regression. 
In each panel are also indicated the average iron content of the 
bin, the root mean square of the linear regression and the number 
of stars. 

The bottom panels of Figures 8 and 9 show the magnitude 
residuals S{M) in the V and K band, respectively, as a function 
of [Fe/H]. These magnitude residuals are calculated as the dif- 
ference between the observed absolute magnitude and the abso- 
lute magnitude as determined from its period using a standard 
PL relation, namely the PLy from Freedman et al. (2001) and 
the the PLk band Persson et al. (2004). In practice, 6{M) is the 
correction to be applied to a universal, [Fe/H] -independent PL 
relation to take into account the effects of metallicity. A positive 
S means that the actual luminosity of a Cepheid is fainter than 
the one obtained with the standard PL relation. The filled squares 
in the bottom panels of Figures 8 and 9 represent the mean val- 
ues of 6(M) in each metallicity bin, the vertical error-bars are 
the errors on the mean and the solid line is the null value, which 
corresponds to an independence of the PL relation on the iron 
content. The horizontal bars indicate the size of the metallicity 
bins. 

Data plotted in the bottom right panel of Fig. 8 indicate 
that the magnitude residuals of the V band in the individual 
bins are positive and located at ^ 2cr (metal-poor) and « 9cr 
(metal-rich) from zero. Moreover, the 6{Mv) in the two bins dif- 
fer from each other at the 3cr level. Our data, then, suggest that 
metal-rich Cepheids in the V band are, at a fixed period, fainter 
than metal-poor ones. Data plotted in the bottom left panel in- 
dicate that this finding is marginally affected by uncertainties 
in the Galactic Cepheid distance scale. The use of Cepheid dis- 
tances provided by Groenewegen (2008) provides a very similar 
trend concerning the metallicity effect. Moreover, the magnitude 
residuals based on the "Old Sample" are located » 1.5cr (metal- 
poor) and X! Act (metal-rich) from zero and the difference in the 
two bins is larger than one cr. 

The results for the K band are displayed in Fig. 9. The 
near-infrared data are on different photometric systems, there- 
fore, we transformed them into the 2MASS (Two Micron All 
Sky Survey) photometric system. In particular, the apparent K 
magnitudes listed in Table 2 and 3 were transformed using the 



Romaniello et al.: The impact of chemical composition on Cepheid properties 



13 







: [Fe/H] = -0.52 / 


/I ; 


: [Fe/H]=0.04 


V 




1 


- rms=0.28 ^, r/ 

- [31 stars] */* 




- rms=0.21 

- [30 stars] ^ ^ 


^ - 


iajO 
6 


1 


• / • 




• jr 
~ • /^ 


- 


> 


1 








- 




CO 

1 


- ,/t 

"/i 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 - 


-/• 1 1 
/I 1 1 1 1 1 1 1 1 1 1 1 1 


. 1 - 



1 1.5 2 1 1.5 2 

log(P) [days] 



0.4 



0.2 



= 0.0 



-0.2 



-0.4 




0.4 


h 


II 


H 


0.2 


- 




- 




1 r 




- 




1 ' 


■ 




J 


L 




- 


0.0 












- 


0.2 


- 




- 


0.4 


New BW distances 







-0.5 

[Fe/H] 



-0.5 

[Fe/H] 



Fig. 8. The top panels show the PL relations calculated in each bin for the V band. The bottom panels show the residuals 6{Mv) as 
a function of the iron abundance for both the "Old" and the "New" Galactic Cepheid distances. The mean values of 6{Mv) in each 
metallicity bin are plotted as filled squares, with the vertical error-bars representing its associated error. The horizontal bars indicate 
the dimension of the bins. The horizontal solid lines indicate the null value which corresponds to an independence of the PL relation 
from the iron content. 



transformations provided by Koen et al. (2007), while the PLk 
relation by Persson was transformed using the transformations 
provided by Carpenter (2001). The correction are typically of 
the order of ~ 0.02 mag. Data plotted in the bottom right panel 
of Fig. 9 indicate that the metal-poor bin is within Icr consis- 
tent with zero. On the other hand, the metal-rich bin differs from 
the null hypothesis by x Act. Moreover, the magnitude residu- 
als in the two metallicity bins differ by x 2cr. This finding is 
at odds with current empirical (Persson et al. 2004; Fouque et 
al. 2007) and theoretical (Bono et al. 1999,2008; Marconi et al. 
2005) evidence suggesting that the PLk relation is marginally 
affected by metal abundance. In order to constrain this effect we 
also adopted the "Old Sample" distances. The magnitude resid- 
uals plotted in the bottom right panel of Fig. 9 show that the two 
metallicity bins are, within the errors, consistent with no metal- 
licity effect. Unfortunately, current error budget (distance, metal 
abundance) does not allow us to reach firm conclusions concern- 
ing the metallicity dependence of the PLk relation. 



5.1. Estimation of the errors caused by sampling of ttie PL 

The error on the residuals showed in Fig. 8 and Fig. 9 is char- 
acteristic of the particular sample of stars we have used. In this 
section we describe the simulations we have performed in order 
to assess the impact of different samplings of the PL relation. 

In order to do so, we have extracted random subsamples 
composed by different numbers of stars from a sample of ob- 
served Cepheid that populate well the PL relation. We choose 
as reference a sample of 771 LMC Cepheids from the OGLE 
database (Udalski et al. 1999) for the test in the V-band and 
the sample of 92 LMC Cepheids (Persson et al. 2004) for the 
test in the /T-band. The latter is the largest observed sample of 
Cepheids in the near-infrared bands. For each extraction a linear 
regression was performed to derive the slope and the zero point 
of the resulting PL relation. Also, we calculated the magnitude 
residuals 6(Mv) and 6{Mk) as defined above for the actual ob- 
served programme stars. We have, then, compared these quanti- 
ties with those derived for the whole sample in order to estimate 
the error due to random sampling and to optimize the number of 
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Fig. 9. The top panels show the PL relations calculated in each bin for the K band. The bottom panels show the residuals 6{Mv) as 
a function of the iron abundance for both the "Old" and the "New" Galactic Cepheid distances. The mean values of 6{Mk) in each 
metallicity bin are plotted as filled squares, with the vertical error-bars representing its associated error The horizontal bars indicate 
the dimension of the bins. The horizontal solid lines indicate the null value which corresponds to an independence of the PL relation 
from the iron content. 



bins we can divide our sample into. It is worth noticing here that 
this procedure will somewhat overestimate the sampling error 
This is because some of the random extractions will generate a 
subsample that covers a small range of periods, while we have 
carefully selected our sample in order to cover a broad period 
range. 

After performing 1 million extractions in each of the bands, 
V and K, for several bin sizes, we have settled for two metallic- 
ity bins of about 30 stars each. With this choice, 6(Mv) results 
to being less that 0.1 mag in 95% of the cases and never larger 
than 0.15. The mean value of the distribution is 0.03 mag. These 
results imply that the non-zero result for the high metallicity bin 
in Figure |8] cannot be due to insufficient sampling of the PL re- 
lation. 

As for the K band, the simulations indicate that S{Mk) is 
smaller that 0.04 mag in 95% of the cases and never larger than 
0.06 mag, with a mean of 0.01 mag. Also in this case, then, the 
results discussed above and displayed in Figure |9] are not signif- 
icantly affected by sampling errors. 



5.2. Comparison with previous results 

We compare our results with two different behaviours (see 
Fig. 10) as examples of the effects of the metallicity on the PL re- 
lation currently available in the literature: independence from the 
iron content and a monotonic decreasing trend (e.g. Kennicutt 
et al. 1998; Sakai et al. 2004; Storm et al. 2004; Groenewegen 
et al. 2004; Macri et al. 2006; Sandage & Tammann 2008), in 
the sense that metal-rich Cepheids are brighter than metal-poor 
ones. 

We use for the comparison the classical results of Kennicutt 
et al. (1998) adopted by the HST Key Project to determine ex- 
tragalactic distances with Cepheids (Freedman et al. 2001). They 
have analysed two Cepheid fields in MlOl, with average values 
of metallicity around -0.4 dex and 0.28 dex (determined from 
measurements of oxygen in H II regions in the two fields). They 
have observed 29 Cepheids in the outer field (low metallicity) 
and 6 1 Cepheids in the inner field (high metallicity) with periods 
between 10 and 60 days. Considering the outcome of our sim- 
ulations, the two Cepheid samples observed by Kennicutt and 
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Table 9. Stellar parameters and Fei, Fen abundances of our Magellanic Cepheids. The Fei values were adopted as final [Fe/H]. 
When available the value from previous studies is also reported ([Fe/H]/,, where L stands for Luck & Lambert 1992; Luck et al. 
1998). In the last column are listed the same values (as in column # 7), but rescaled to the iron solar abundance adopted in our work. 



ID 


Teff 


V, 


logg 


[Fei/H] 


[Feii/H] 


[Fe/H], 


[Fe/H]i, 


LMC 


HV877 


4690 


5.40 


0.5 


-0.44 


-0.47 






HV879 


5630 


3.05 


1.0 


-0.14 


-0.14 


-0.56* 


-0.46 


HV971 


5930 


2.30 


1.4 


-0.29 


-0.29 






HV997 


5760 


3.10 


1.2 


-0.21 


-0.22 






HV 1013 


4740 


5.35 


0.2 


-0.59 


-0.60 






HV 1023 


5830 


3.10 


1.1 


-0.28 


-0.27 






HV 2260 


5770 


3.40 


1.6 


-0.38 


-0.36 






HV 2294 


5080 


3.90 


0.5 


-0.42 


-0.42 


-0.06° 


+0.10 


HV 2337 


5560 


3.30 


1.6 


-0.35 


-0.36 






HV 2352 


6100 


3.65 


1.6 


-0.49 


-0.47 






HV 2369 


4750 


6.00 


0.3 


-0.62 


-0.63 






HV 2405 


6170 


4.20 


2.3 


-0.27 


-0.28 






HV 2580 


5360 


2.75 


0.7 


-0.24 


-0.25 






HV 2733 


5470 


2.90 


1.8 


-0.28 


-0.27 






HV 2793 


5430 


2.90 


0.9 


-0.10 


-0.11 






HV 2827 


4790 


4.00 


0.0 


-0.38 


-0.33 


-0.24* 


-0.14 


HV 2836 


5450 


2.85 


1.0 


-0.16 


-0.19 






HV 2864 


5830 


2.80 


1.5 


-0.19 


-0.20 






HV 5497 


5100 


3.40 


0.3 


-0.25 


-0.24 


-0.48* 


-0.38 


HV 6093 


5790 


4.50 


1.5 


-0.60 


-0.60 






HV 12452 


5460 


2.90 


1.0 


-0.35 


-0.37 






HV 12700 


5420 


3.15 


1.4 


-0.36 


-0.35 






SMC 


HV817 


5850 


3.25 


1.0 


-0.82 


-0.84 






HV823 


5990 


3.80 


1.4 


-0.80 


-0.81 






HV824 


5170 


3.00 


0.7 


-0.73 


-0.74 


-0.94* 


-0.84 


HV829 


5060 


3.30 


0.2 


-0.76 


-0.73 


-0.61* 


-0.51 


HV834 


5750 


2.95 


1.2 


-0.63 


-0.64 


-0.59* 


-0.49 


HV837 


5140 


2.90 


0.0 


-0.83 


-0.80 


-0.91* 


-0.81 


HV847 


4790 


2.80 


0.0 


-0.75 


-0.77 






HV865 


6130 


1.90 


0.5 


-0.87 


-0.82 


-0.44° 


-0.28 


HV 1365 


5340 


2.48 


0.6 


-0.82 


-0.84 






HV 1954 


5890 


2.47 


1.0 


-0.76 


-0.75 






HV 2064 


5550 


3.10 


0.7 


-0.64 


-0.64 






HV 2195 


5970 


2.90 


1.0 


-0.67 


-0.68 


-0.45° 


-0.29 


HV 2209 


6130 


2.30 


1.1 


-0.65 


-0.67 






HV 11211 


4830 


2.60 


0.0 


-0.83 


-0.81 






" Luck & Lambert (1992). 












* Luck et al. 


(1998). 















collaborators are comparable to our sub-samples in each bin, it is 
then reasonable to compare the qualitative indications about the 
effect of the metallicity that can be derived from the two analy- 
ses. The complete comparison is possible only in the V-band. 

As we have already mentioned the error due to the sampling 
of the PL relation is much smaller than the one on the residuals. 
Data plotted in Figures 8,9 and 10 disclose several circumstan- 
tial evidence. This means that the increasing trend of the latter 
in the V-band, as a function of the iron content, is real with a 
confidence level of 99%. 

i) V-band case 

- No dependence of S{Mv) on [Fe/H]: A null effect on the 
metal abundance would imply that the residuals of the 
two bins should be located, within the errors, along the 
independence line. Current findings and the outcome of the 
simulations mentioned above indicate that this hypothesis 
can be excluded completely. The increasing trend of the 



residual in the V-band, as a function of the iron content, is 
real with a confidence level larger than 95%. 



Monotonically decreasing diMy): we compare the classical 
results of Kennicutt et al. (1998, open circles and solid 
line in Fig. 10) with our data. Since the increasing trend of 
our residuals is real, this hypothesis is incompatible with 
our results. In passing, we also note that the reddening 
estimates adopted by Kennicutt et al. (1998) are based on 
the observed mean colors. However, metal-poor Cepheids 
are, at fixed period, hotter than metal-rich ones. Therefore, 
Kennicutt et al. would have underestimated the reddenings 
and the luminosities of his metal-poor sample, producing an 
apparent underluminosity for metal-poor Cepheids. 



//) /T-band case 
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Fig. 10. The V-band residuals compared to Freedman et al. 
(2001) PL relation are plotted against the iron content measured 
from observed spectra (bottom right panel in Fig. 8). The filled 
squares display the mean value in each metallicity bin, with 
its associated errorbar The metallicity dependence estimated by 
Kennicutt et al. (1998) using two Cepheid fields in MlOl (open 
circles) is shown as a full line. 



- Current data do not allow us to reach a firm conclusion 
concerning the metallicity effect. 



To summarize, we found an increasing trend of the V-band 
residuals with the iron content. This result is in disagreement 
with an independence of the PL relation on iron abundance and 
with the linearly decreasing trend found by other observational 
studies in the literature (e.g. Kennicutt et al. 1998). 

5.3. LMC distance: "short" scale vs "long" scale 

Regrettably we can find in the recent literature values of the 
LMC distance modulus ranging from 18.1 to 18.8, not always 
obtained with different techniques. Those studies finding a dis- 
tance modulus less than 18.5 support the so-called "short" dis- 
tance scale, whereas those finding it greater than 18.5 support 
the "long" distance scale (for a review of the results and methods 
see Benedict et al. 2002; Gibson 2000). More recently, Schaefer 
(2008) found that distance estimates to the LMC published be- 
fore 2001 present a large spread ( 18.1 < i-i < 18.8). On the other 
hand, distances published after 2002 tightly concentrate around 
the value adopted by the HST Key Projects (/i - 18.5, Freedman 
et al. 2001; Saha et al. 2006). In order to overcome this suspi- 
cious bias, we decided to investigate the different behaviour of 
the PL relation depending on the adopted LMC distance scale. In 
order to do that, we have repeated our calculation of the V-band 
residuals assuming a distance modulus for LMC of 18.3 (repre- 
sentative of the "short" scale) and of 18.7 (for the "long" scale). 
The results are shown in Fig. 1 1 . 

In the left and right panel, the V-band residuals refer to the 
"short" and the "long" distance scale, respectively. For /ilmc of 
18.3, the V-band residuals are located at least at 5 cr from zero. 
The difference between the metal-poor and the metal-rich bin is 



of the order of one cr. This trend is similar to the trend we ob- 
tained using yU£MC=18.5, however, the distance from zero of the 
metal-poor bin is significantly larger. This result disagrees with 
an independence of the PL relation on the iron abundance and the 
monotonic decreasing behaviour at the 99.99 % level (according 
to the;t'^ method). On the other hand, the 6(Mv) values for fiiMC 
of 18.7 are located at ^ 2 (metal-poor) and ~ 8 (metal-rich) cr 
from zero. The difference between the two bins is at least at 4 
cr level. The data trend is slightly steeper than for //£mc = 18.5. 
Using again a x^ technique, we find that this result disagrees 
with an independence of the PL relation on iron abundance and 
with the linearly decreasing trend often quoted in the literature 
(e.g. Kennicutt et al. 1998). 

The results based on the tests performed assuming different 
LMC distances are the following: 

Short scale - Data plotted in the left panel of Fig. 1 1 indicate that 
V-band PL relation does depend on the metal content. Indeed, 
the two bins are located at least at 5 cr from zero. This means 
that the zero-point of the quoted PL relations do depend on the 
metal abundance. However, the difference between the two bins 
is small and of the order of one cr. This indicates that the metal- 
licity effect supported by the short scale is mainly caused by a 
difference in the zero-point. This result would imply a significant 
difference in the zero-point of Magellanic Cepheids. However, 
such a difference is not supported by current empirical (Laney 
& Stobie 2004; van Leeuwen et al. 2007; Fouque et al. 2007; 
Sandage & Tammann 2008) and theoretical evidence. 

Long scale - Data plotted in the right panel of Fig. 1 1 indi- 
cates that V-band PL relation does depend on the metal content. 
The difference from zero ranges from ~ 2 for the metal-poor bin 
to more than 8 cr for the metal-rich bin. Moreover, the two bins 
differ at 4.5cr level. This means that the metallicity effect sup- 
ported by the long scale is mainly caused by a difference in the 
slope. This finding, taken at face value, would imply that metal- 
poor Cepheids (e.g. Cepheids in IC 1613, -1.3 < [Fe/H] < 
-0.7, Skillman et al. 2003) should be, at fixed period, « 0.5 (V) 
mag brighter than Galactic Cepheids. This difference is not sup- 
ported by current empirical (Dolphin et al. 2003; Antonello et 
al. 2006; Pietrzynski et al. 2006; Saha et al. 2006; Fouque et al. 
2007) and theoretical evidence. 

The quoted results suggest that the V-band PL relation is af- 
fected by metal abundance. This finding is marginally affected 
by the adopted LMC distance. It is worth mentioning that recent 
empirical estimates based on robust primary standard candles 
indicate that the true LMC distance is 18.5 + 0.1 (Alves 2004; 
Benedict el al. 2007; van Leeuwen et al. 2007; Catelan & Cortes 
2008; Feast et al. 2008; Groenewegen et al. 2008; Sollima et al. 
2008). In view of this convergence on the LMC distance and on 
the results based on the "short" and on the "long" distance scale, 
the results based on fiiMC = 18.5 appear to be the most reliable 
ones. 

6. Conclusions 

We have directly measured the iron abundances for 68 Galactic 
and Magellanic Cepheids from FEROS and UVES high resolu- 
tion and high signal-to-noise spectra. We have used these mea- 
surements to assess the influence of the stellar iron content on the 
Cepheid PL relation in the V and in the K band. In order to do 
this we have related the V-band and the TT-band residuals from 
the standard PL relations of Freedman et al. (2001) and Persson 
et al. (2004), respectively, to [Fe/H]. Differently from previ- 
ous studies, we can constrain the PL relation using Cepheids 
with known distance moduli and chemical abundances, homo- 
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Fig. 11. The V-band residuals compared to the Freedman et al. (2001) PL relation versus the iron content measured from observed 
spectra, assuming an LMC distance modulus of 18.3 (left) and of 18.7 (right). The filled squares represent the mean value in each 
metallicity bin, with its associated errorbar. 



geneously measured, that cover almost a factor of ten in metal- 
licity. 

For our Galactic sample, we find that the mean value of the 
iron content is solar (cr = 0.10, see Fig. 6), with a range of values 
between -0.18 dex and -1-0.25 dex. For the LMC sample, we find 
that the mean value is about ~ -0.33 dex (cr =0.13, see Fig. 6), 
with a range of values between -0.62 dex and -0. 10 dex. For the 
SMC sample, we find that the mean value is about ~ -0.75 dex 
(o- - 0.08, see Fig. 6), with a range of values between -0.87 and 
-0.63. 

We have compared our results with the analyses of FC97, 
Andrievsky et al. (2002a, 2002b, 2002c) and Luck et al. 
(2003) for the Galactic Cepheids and LL92 and L98 for the 
Magellanic Clouds. Regarding the Galactic sample, our results 
are marginally more in agreement with Andrievsky's values than 
with FC97 and the diff'erences, on average, appear rather small. 
Considering the Magellanic Cepheids, we have a poor agreement 
with LL92, which could be in part accounted for by different an- 
alytical tools and data quality. Our data are in better agreement 
with L98 results and the spread of our iron abundances in the 
LMC and SMC is similar to the one they reported. We note that 
the mean metallicity that they found with their complete sample 
(-0.30 dex and -0.74 dex) is in very good agreement with our 
results. 

Our main results concerning the effect of the iron abundance 
on the PL relation are summarized in Fig. 8 and Fig. 9 (bottom 
panels) and they hold for a LMC distance modulus of 18.50. In 
Fig. 10 is also showed the comparison in the V-band with the 
empirical results of Kennicutt et al. (1998) in two Cepheid fields 
of MlOl (open circles and solid line). The main findings can be 
summarized as follows: 



- The V-band PL relation does depend on the metal abun- 
dance. This finding is marginally affected by the adopted 
distance scale for the Galactic Cepheids and by the LMC 
distance. 



- Current data do not allow us to reach a firm conclusion 
concerning the dependence of the /T-band PL relation on 
the metal content. The use of the most recent distances for 
Galactic Cepheids (Benedict et al. 2007; Fouque et al. 2007; 
van Leeuwen et al. 2007) indicates a mild metallicity effect. 
On the other hand, the use of the old distances (Storm et al. 
2004) suggest a vanishing effect. 

- Residuals based on a canonical LMC distance (;U/,MC = 18.5) 
and on the most recent distances for Galactic Cepheids 
present a well defined effect in the V-band. The metal-poor 
and the metal-rich bin are x; 2cr and » 9cr from the null hy- 
pothesis. Moreover, the two metallicity bins differ at the 3cr 
level. 

- By assuming a "short" LMC distance (fJ-iMC - 18.3) the 
residuals present a strong metallicity dependency in the zero- 
point of the V-band PL relation. By assuming a "long" LMC 
distance (jilmc = 18.7) we found a strong metallicity effect 
when moving from metal-poor to metal-rich Cepheids. This 
indicates a significant change in the slope and probably in 
the zero-point. The findings based on the "short" and on the 
"long" LMC distance are at odds with current empirical and 
theoretical evidence, suggesting a smaller metallicity effect. 

- Metal-rich Cepheids in the V-band are systematically fainter 
than metal-poor ones. This evidence is strongly supported by 
the canonical, the "short" and the "long" LMC distance. 

The above results together with recent robust LMC distance 
estimates indicate that the behaviours based on the canonical dis- 
tance appear to be the most reliable ones. 

In order to constrain on a more quantitative basis the metal- 
licity dependence of both zero-point and slope of the optical 
PL relations is required a larger number of Cepheids covering 
a broader range in metal abundances. Moreover, for each metal- 
licity bin Cepheids covering a broad period range are required to 
reduce the error on the residuals and to constrain on a quantita- 
tive basis the fine structure of the PL relation in optical and NIR 
photometric bands. 
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Thesaurus of key words used in the annual subject 
indexes (valid from January 1997) 

The list is common to Astronomy and Astrophysics, The 
Astrophysical Journal and Monthly Notices of the Royal 
Astronomical Society. In order to ease the search, the key 
words are subdivided into broad categories. 

No more than six codes all together should be listed for 
a paper as this is the limit fixed by the computer program. 

The key words in boldface listed under the code num- 
bers 07.03.2, 07.16.2, 08.09.2, 08.16.5, 08.19.5, 09.09.1, 
09.16.2, 10.07.3, 10.15.2, 11.02.2, 11.09.1, 11.17.4 are in- 
tended for use with specific astronomical objects; each con- 
tains the word "individual" . The corresponding code num- 
bers should never be used alone, but always in combination 
with the most common names for the astronomical objects 
in question. For example, if a paper discusses three individ- 
ual galaxies, these should be coded on the title page of the 
manuscript in the following manner: 

11.09.1 Arp 220; 11.09.1 M 51; 11.09.1 NGC 4472 

Note that each object (in the example the three galaxies) 
counts as one code within the allowed limit of six codes. 

The parts of the key words in italics are for reference 
only and should be omitted when the key words are entered 
on the manuscript. 

General 

01.02.1 Book reviews 

01.05.1 Editorials notices 

01.05.2 Errata, addenda 

01.05.3 Extraterrestrial intelligence 
01.08.1 History and philosophy of astronomy 
01.13.1 Miscellaneous 

01.15.1 Obituaries, biographies 

Physical data and processes 

02.01.1 Acceleration of particles 

02.01.2 Accretion, accretion disks 

02.01.3 Atomic data 

02.01.4 Atomic processes 
02.02.1 Black hole physics 

02.03.1 Chaos 

02.03.2 Conduction 

02.03.3 Convection 

02.03.4 Cosmic strings 

02.04.1 Dense matter 

02.04.2 Diffusion 

02.05.1 Elementary particles 

02.05.2 Equation of state 

02.07.1 Gravitation 

02.07.2 Gravitational waves 
02.08.1 Hydrodynamics 
02.09.1 Instabilities 

02.12.1 Line: formation 

02.12.2 Line: identification 

02.12.3 Line: profiles 

02.13.1 Magnetic fields 

02.13.2 Magnetohydrodynamics (MHD) 

02.13.3 Masers 

02.13.4 Molecular data 

02.13.5 Molecular processes 



02.14.1 Nuclear reactions, nucleosynthesis, abundances 

02.16.1 Plasmas 

02.16.2 Polarization 

02.18.5 Radiation mechanisms: non-thermal 

02.18.6 Radiation mechanisms: thermal 

02.18.7 Radiative transfer 

02.18.8 Relativity 
02.19.2 Scattering 
02.19.1 Shockwaves 
02.20.1 Turbulence 

02.23.1 Waves 

Astronomical instrumentation, methods and tech- 
niques 

03.01.2 Atmospheric effects 
03.02.1 Balloons 

03.09.1 Instrumentation: detectors 

03.09.2 Instrumentation: interferometers 

03.09.3 Instrumentation: miscellaneous 

03.09.4 Instrumentation: photometers 

03.09.5 Instrumentation: polarimeters 

03.09.6 Instrumentation: spectrographs 

03.13.1 Methods: analytical 

03.13.2 Methods: data analysis 

03.13.7 Methods: laboratory 

03.13.3 Methods: miscellaneous 

03.13.4 Methods: numerical 

03.13.5 Methods: observational 

03.13.6 Methods: statistical 

03.19.1 Site testing 

03.19.2 Space vehicles 

03.20.1 Techniques: image processing 

03.20.2 Techniques: interferometric 

03.20.3 Techniques: miscellaneous 

03.20.4 Techniques: photometric 

03.20.5 Techniques: polarimetric 

03.20.6 Techniques: radar astronomy 

03.20.7 Techniques: radial velocities 

03.20.8 Techniques: spectroscopic 

03.20.9 Telescopes 

Astronomical data bases 



04.01.1 


Astronomical data bases: miscellaneous 


04.01.2 


Atlases 


04.03.1 


Catalogs 


04.19.1 


Surveys 



Astrometry and celestial mechanics 



05.01.1 


Astrometry 


05.03.1 


Celestial mechanics, stellar dynamics 


05.05.1 


Echpses 


05.05.2 


Ephemerides 


05.15.1 


Occultations 


05.18.1 


Reference systems 


05.20.1 


Time 



The Sun 




06.01.1 


Sun: abundances 


06.01.2 


Sun 


activity 


06.01.3 


Sun 


atmosphere 


06.03.1 


Sun 


chromosphere 


06.03.2 


Sun 


corona 


06.05.1 


Sun 


evolution 


06.06.1 


Sun 


faculae, plages 


06.06.2 


Sun 


filaments 


06.06.3 


Sun 


flares 


06.06.4 


Sun 


fundamental parameters 


06.07.1 


Sun 


general 


06.07.2 


Sun 


granulation 


06.09.2 


Sun 


infrared 


06.09.1 


Sun 


interior 


06.13.1 


Sun 


magnetic fields 


06.15.1 


Sun 


oscillations 


06.16.1 


Sun 


particle emission 


06.16.2 


Sun 


photosphere 


06.16.3 


Sun 


prominences 


06.18.1 


Sun 


radio radiation 


06.18.2 


Sun 


rotation 


06.19.1 


(Sun:) solar-terrestrial relations 


06.19.2 


(Sun:) solar wind 


06.19.3 


(Sun:) sunspots 


06.20.1 


Sun: transition region 


06.21.1 


Sun: UV radiation 


06.24.1 


Sun: X-rays, gamma rays 


Solar system 


07.03.1 


Comets: general 


07.03.2 


Comets: individual: . . . 


07.05.1 


Earth 


07.09.1 


Inteplanetary medium 


07.13.1 


Meteors, meteoroids 


07.13.2 


Minor planets, asteroids 


07.13.3 


Moon 


07.16.1 


Planets and satellites: general 


07.16.2 


Planets and satellites: 




individual: . . . (alphabetic order) 


07.19.1 


Solar system: formation 


07.19.2 


Solar system: general 


Stars 




08.01.1 


Stars: abundances 


08.01.2 


Stars: activity 


08.16.4 


Stars: AGB and post-AGB 


08.01.3 


Stars: atmospheres 


08.02.1 


(Stars:) binaries (including multiple): close 


08.02.2 


(Stars:) binaries: eclipsing 


08.02.3 


(Stars:) binaries: general 


08.02.4 


(Stars:) binaries: spectroscopic 


08.02.5 


(Stars:) binaries: symbiotic 


08.02.6 


(Stars:) binaries: visual 


08.02.7 


(Stars:) blue stragglers 


08.03.1 


Stars: carbon 


08.03.2 


Stars: chemically pecuhar 


08.03.3 


Stars: chromospheres 


08.03.4 


(Stars:) circumstellar matter 


08.03.5 


Stars: coronae 


08.04.1 


Stars: distances 


08.05.1 


Star 


s: early-type 



08.05.2 
08.05.3 
08.06.1 
08.06.2 
08.06.3 



08.07.1 
08.08.1 
08.08.2 
08.09.1 
08.09.2 
08.09.3 
08.11.1 
08.12.1 
08.12.2 
08.12.3 
08.13.1 
08.13.2 
08.14.1 
08.14.2 
08.15.1 
08.16.1 
08.16.2 
08.16.3 
08.16.5 
08.16.6 
08.16.7 
08.18.1 
08.19.1 
08.19.2 
08.19.3 
08.19.4 
08.19.5 
08.19.6 
08.22.1 
08.22.2 
08.22.3 
08.23.1 
08.23.2 



Stars: emission-line. Be 

Stars: evolution 

Stars: flare 

Stars: formation 

Stars: fundamental parameters 

(classification, colors, luminosities, masses, 

radii, temperatures, etc.) 

Stars: general 

(Stars:) Hertzsprung-Russel (HR) diagram 

Stars: horizontal-branch 

Stars: imaging 

Stars: individual: . . . 

Stars: interiors 

Stars: kinematics 

Stars: late-type 

Stars: low-mass, brown dwarfs 

Stars: luminosity function, mass function 

Stars: magnetic fields 

Stars: mass-loss 

Stars: neutron 

(Stars:) novae, cataclysmic variables 

Stars: oscillations (including pulsations) 

Stars: peculiar (except chemically peculiar) 

(Stars:) planetary systems 

Stars: Population 11 

Stars: pre-main sequence 

(Stars:) pulsars: general 

(Stars:) pulsars: individual: . . . 

Stars: rotation 

Stars: statistics 

(Stars:) subdwarfs 

(Stars:) supergiants 

(Stars:) supernovae: general 

(Stars:) supernovae: individual: . . . 

(Stars:) starspots 

(Stars: variables:) Cepheids 

(Stars: variables:) 6 Scu 

Stars: variables: other 

(Stars:) white dwarfs 

Stars: Wolf-Rayet 



Interstellar medium (ISM), nebulae 

09.01.1 ISM: abundances 

09.01.2 ISM: atoms 
09.02.1 ISM: bubbles 

09.03.1 ISM: clouds 

09.03.2 (ISM:) cosmic rays 
09.04.1 (ISM:) dust, extinction 
09.07.1 ISM: general 

09.08. 1 (ISM:) H ii regions 

09.09.1 ISM: individual objects: . . . 

(except planetary nebulae) 

09.10.1 ISM: jets and outflows 

09.11.1 ISM: kinematics and dynamics 

09.13.1 ISM: magnetic fields 

09.13.2 ISM: molecules 

09.16.1 (ISM:) planetary nebulae: general 

09.16.2 (ISM:) planetary nebulae: individual: 
09.18.1 ("/S'M.-) reflection nebulae 

09.19.1 ISM: structure 

09.19.2 ISM: supernova remnants 



The Galaxy 

10.01.1 Galaxy: abundances 

10.03.1 Galaxy: center 

10.05.1 Galaxy: evolution 

10.06.1 Galaxy: formation 

10.06.2 Galaxy: fundamental parameters 

10.07.1 Galaxy: general 

10.07.2 (Galaxy:) globular clusters: general 

10.07.3 (Galaxy:) globular clusters: individual: . . . 
10.08.1 Galaxy: halo 

10.11.1 Galaxy: kinematics and dynamics 

10.15.1 (Galaxy:) open clusters and associations: general 

10.15.2 (Galaxy:) open clusters and associations: in- 
dividual: . . . 

10.19.1 (Galaxy:) solar neighbourhood 

10.19.2 Galaxy: stellar content 

10.19.3 Galaxy: structure 

Galaxies 

11.01.1 Galaxies: abundances 

11.01.2 Galaxies: active 

11.02.1 (Galaxies:) BL Lacertae objects: general 

11.02.2 (Galaxies:) BL Lacertae objects: individual: 

11.03.1 Galaxies: clusters: general 

11.03.4 Galaxies: clusters: individual: ... 

11.03.2 Galaxies: compact 

11.03.3 (Galaxies:) cooling flows 
11.04.1 Galaxies: distances and redshifts 

11.05.1 Galaxies: elliptical and lenticular, cD 

11.05.2 Galaxies: evolution 

11.06.1 Galaxies: formation 

11.06.2 Galaxies: fundamental parameters 
(classification, colors, luminosities, masses, 
radii, etc.) 

11.07.1 Galaxies: general 

11.08.1 Galaxies: halos 

11.09.1 Galaxies: individual: ... 

11.09.2 Galaxies: interactions 

11.09.3 (Galaxies:) intergalactic medium 

11.09.4 Galaxies: ISM 

11.09.5 Galaxies: irregular 
11.10.1 Galaxies: jets 

11.11.1 Galaxies: kinematics and dynamics 

11.12.1 (Galaxies:) Local Group 

11.12.2 Galaxies: luminosity function, mass function 

11.13.1 (Galaxies:) Magellanic Clouds 

11.13.2 Galaxies: magnetic fields 
11.14.1 Galaxies: nuclei 

11.16.1 Galaxies: photometry 

11.16.2 Galaxies: peculiar 

11.17.1 (Galaxies:) quasars: absorption lines 

11.17.2 (Galaxies:) quasars: emission lines 

11.17.3 (Galaxies:) quasars: general 

11.17.4 (Galaxies:) quasars: individual: ... 

11.19.1 Galaxies: Seyfert 

11.19.2 Galaxies: spiral 

11.19.3 Galaxies: starbust 

11.19.4 Galaxies: star clusters 
11.19.7 Galaxies: statistics 

11.19.5 Galaxies: stellar content 

11.19.6 Galaxies: structure 



Cosmology 

12.03.1 (Cosmology:) cosmic microwave background 

12.03.2 Cosmology: miscellaneous 

12.03.3 Cosmology: observations 

12.03.4 Cosmology: theory 

12.04.1 (Cosmology:) dark matter 

12.04.2 (Cosmology:) diffuse radiation 

12.04.3 (Cosmology:) distance scale 
12.05.1 (Cosmology:) early Universe 
12.07.1 (Cosmology:) gravitational lensing 

12.12.1 (Cosmology:) large-scale structure of Universe 

Sources as a function of wavelength 

13.07.1 Gamma rays: bursts 

13.07.2 Gamma rays: observations 

13.07.3 Gamma rays: theory 

13.09.1 Infrared: galaxies 

13.09.2 Infrared: general 

13.09.3 Infrared: ISM: continuum 

13.09.4 Infrared: ISM: lines and bands 

13.09.5 Infrared: solar system 

13.09.6 Infrared: stars 

13.18.1 Radio continuum: galaxies 

13.18.2 Radio continuum: general 

13.18.3 Radio continuum: ISM 

13.18.4 Radio continuum: solar system 

13.18.5 Radio continuum: stars 

13.19.1 Radio lines: galaxies 

13.19.2 Radio lines: general 

13.19.3 Radio lines: ISM 

13.19.4 Radio lines: solar system 

13.19.5 Radio lines: stars 

13.21.1 Ultraviolet: galaxies 

13.21.2 Ultraviolet: general 

13.21.3 Ultraviolet: ISM 

13.21.4 Ultraviolet: solar system 

13.21.5 Ultraviolet: stars 

13.25.1 X-rays: bursts 

13.25.2 X-rays: galaxies 

13.25.3 X-rays: general 

13.25.4 X-rays: ISM 

13.25.5 X-rays: stars 



